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ULTRASONIC DELAY LINES. I.* 
BY 
H. B. HUNTINGTON,} A. G. EMSLIE t AND V. W. HUGHES,§ 
Massachusetts Institute of Technology, 


Cambridge, Massachusetts. 


ABSTRACT. 


This paper, the first of two dealing with ultrasonic delay lines, is concerned mainly with 
the theoretical considerations involved in the general process of ultrasonic propagation. The 
transducer action is given considerable attention and an equivalent circuit employing a trans- 
mission line as an element is demonstrated. The diffraction spreading of the beam and ab- 
sorptive losses in the transmitting medium are also discussed. 


A. INTRODUCTION. 


This paper and its sequel give an account of that research on ultra- 
sonics which was carried out at the Radiation Laboratory at M. I. T. 
on the development of ultrasonic delay lines. In this paper will be 
presented some of the general theoretical considerations which are 
involved in such a device and in the second paper the construction and 
operation of the actual delay line is treated from the experimental 
standpoint. 

An ultrasonic delay line functions by transforming an electrical sig- 
nal into a sonic one by an electromechanical transducer, propagating it 
through a fixed path in the transmitting medium and then converting 
the sonic signal back into an electrical signal by a second transducer. 
Very long delays can be obtained by this method, in comparison to 
purely electrical devices, because of the greatly decreased velocity of 
sound as compared to the velocity of electromagnetic propagation. 
Since the velocity of sound in liquids is about a third of that in solids, 


* This paper is based on work done for the Office of Scientific Research and Development 
under Contract OEMsr-262 with the Massachusetts Institute of Technology. 

t Now at Rensselaer Polytechnic Institute, Troy, New York. 

t Now at Williams College, Williamstown, Massachusetts. 


§ Now at Columbia University, New York, New York. 
(Note—The Franklin Institute is not responsible for the statements and opinions advanced by contributors in 
the JouRNAL.) 
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a good deal of the research has been on liquid delay lines. However, 
some promising results have been obtained with multiply reflected beams 
in solids. 

The frequency response of an ultrasonic delay line has a maximum 
at{the resonant frequency of the quartz crystals and falls off at low fre- 
quencies. Therefore, it is necessary, in order to obtain undistorted 
transmission of pulses, to send the pulses through the line on a carrier 
whose frequency is near that of the crystals. This allows exceedingly 
faithful transmission. In some cases, however, the wave shape of the 
transmitted pulse is not so critical provided a rapid rate of rise is main- 
tained. This is true when it is desired to delay a trigger pulse and there- 
fore a trigger line can be operated without a carrier. 


B. TRANSDUCER ACTION AND LOSS THROUGH MISMATCH. 


The treatment of this section deals with the analysis of the funda- 
mental action of the piezoelectric transducer, operating in the thickness 
vibration. The results deduced here are subsequently referred in 
another paper and form the basis for approach to many design problems 
in ultrasonic devices. As such these results are by no means original 
with the present group but the presentation differs somewhat from those 
already in the literature. The general plan has been to attack the 
problem directly by fundamental procedure, making use of assumptions 
and approximations valid for the situations under consideration but 
avoiding equivalent circuits which hold only over a limited range of 
frequency. From fundamental equations an expression for the im- 
pedance of the transducer is obtained as seen from the external driving 
circuit. Transformation of the impedance in series with the electro- 
static capacity to a shunt admittance is then carried out and certain 
small terms are neglected. The behavior of the circuit at resonance is 
investigated. 

At this point a somewhat more general treatment is introduced to 
show that, if the fundamental equations are written so that they ex- 
press the current-voltage relations in a six-terminal network and the 
same approximation is made on transforming from a series to a shunt 
circuit, there exists a very simple equivalent circuit involving a trans- 
mission line section. This circuit is then a general statement of all 
the specific results which have been derived earlier. 

The next paragraphs investigate the frequency response of the trans- 
ducer. The behavior near resonance is analyzed in terms of an equiva- 
lent “Q.” 

The last part of this section is concerned with the calculation of the 
loss in energy transfer arising from impedance mismatch. For broad 
band performance the electrical circuits on both driving and receiving 
crystals must be loaded with resistors of a few hundred ohms. On the 
other hand, the transducers usually present impedances of the order of 
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several kilohms. Some consideration is given to the choice of crystal 
diameter for least loss, keeping bandwidth constant. The possibilities 
in driving the crystal at higher harmonics are also discussed. 

The first step will be the development of an expression for the elec- 
trical impedance of the transducer in terms of acoustic impedances. 
Consider a piezoelectric crystal of uniform thickness fitted across a 
tube of uniform cross section and bounded by media 1 and 2 respectively, 
at opposite faces (see Fig. 1). 

The case considered is that of vibration propagated in the thickness 
(x) direction as for example, in X-cut quartz. It will be assumed that 
no other modes are appreciably excited and that the displacements de- 
pend only on the x-coordinate. The dimensions of the vibrating crystal 
surfaces are large compared to the wavelength of the vibration. No 
acoustic reflection is returned from objects in the media during the time of 
excitation. One would like to know what is the equivalent electrical 
circuit as seen looking at terminals attached respectively to the elec- 


trodes on the crystal faces. 


e -jAx Pe 
Ce ips 


Medium |! Crystal Medium 2 
Fic. 1. Vibrating crystal loaded on both sides. 


Ae ak B 


» % aa. 
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The appropriate equations for the piezoelectric medium are ! 


dg 
ax 
ae 
eu ae & 


where ~ = displacement of the mechanically vibrating medium, 
Ss = stress, 
Y = the appropriate modulus of elasticity, 
D = the electrical displacement = 47, 
o = the free charge per unit area of the electrodes, 
E = electric field. 
K = dielectric constant of the crystal, 
h = piezoelectric coefficient appropriate to the independent vari- 
‘ ables used and the orientation of the crystal. 


It should be remarked that Y, f and K are functions of the orientation 
of the crystal. 


1 “Electromechanical Transducers and Wave Filters,” W. P. Mason, Van Nostrand, pp. 
202, 231 (1942); NDRC Report 14-271, W. Swann. 
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When the crystal is electrically excited at a frequency w/2z, nearly 
plane waves are propagated outward from both interfaces with wave- 
lengths of 27/8, and 27/82, respectively (see Fig. 1). In the adjacent 
media the stress is given by S = Yi .,2 (d&/dx). Inside the crystal, 
waves will be moving in both directions. (It will be our practice 
throughout these articles to consider that the time dependent factor, 
e/*, is implicitly contained in the amplitude coefficient.) One is able 
to assign a definite wavelength to the propagation inside the crystal by 
virtue of the fact that the quantity hD/4m in (1) is independent of x. 

The boundary conditions are that the mechanical displacement and 
pressure be continuous. This gives four equations for the four & un- 
known amplitude coefficients: 


A=B+C, (3a) 
Be~ #64 +. Ce ibd = A’e~ibrd, (3b) 
jB:Y:A = — j8VYB + jBYC + hD/4r, (3c) 
—jBYBe~ #4 + jBYCei#4 + hD/4x = — jB.V2A'e~ i624, = (3d) 


where (1) is used to obtain the last two equations. 
Taking S as the area of the interface, one writes for the current 


(4) 


where (c) is free change per unit area of the electrodes. For the voltage 
across the crystal measured in the sense of Ohm’s Law, one obtains from 
equation (2) 


= + , Edx = h(éa Bs fo) + sh (5) 


The equations (3) can be solved for displacements in the usual way and 


BY {cos (6d) — 1} + 7Y:6: sin Bd (6a) 
ho YAY, + B:Y2) cos Bd + gL (BY)? + Bi YiB2Y2 | sin Bd 


= jhoF,(8d), 
cae eee | 


BY {cos (Bd) — 1} + 7 Y2f2 sin pd (6b) 
* (BY)(8:Y: + B2Y2) cos Bd + JL(BY)? + B:Y:62¥2] sin Bd 


= jhoF;(8d). 


The expression for the electrical impedance of the crystal can now be 


where 


is the 
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written 
(6c) 


where 
(6d) 


The first term in (6c) is always much larger than the second, pro- 
vided that at least one of the media adjoining the crystal is either a 


() Le 


Fic. 2. Shunt circuit. | 


liquid or a solid. This makes it feasible to express the transducer as a 
simple parallel combination. 


I/V = jwC + (Z)7 (7) 
Here C is the electrostatic capacity of the crystal and 
162°d? 


2 (7a) 


~ RK %wS F(Bd) 


is the equivalent electrical impedance of the crystal under the above 
conditions of acoustical loading (see Fig. 2). The approximation im- 
plicit in (7) greatly simplifies the subsequent treatment and the results 
check closely with the more exact calculation. For example, the 


0 


0 


Fic. 3. Circuit with generator. 


resonant frequency of the transducer as developed here differs by less 
than 1 per cent. from the exact value. 

From (6c) it is apparent that the electrostatic capacity of the crystal 
acts as a series element. This is a consequence of the fact that the 
thickness vibration is being considered. (Where the direction of the 
supersonic waves is at right angles to the electric field the electrostatic 
capacity appears as a shunt element.) 
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If one had developed the case of the receiving rather than the trans- 
mitting crystal, an impinging acoustical wave would have been in- 
troduced. Then the equivalent electrical circuit would have involved 
a generator, with voltage amplitude proportional to the amplitude of 


the incoming sonic wave (see Fig. 3). 
The resonant condition for the crystal is that its width be an odd 


number of half wavelengths, or Bd = mx where n is odd. At these 
values F (see 6) reduces to 4(8:Y: + 62Y:)—. Since Y = v*p where p 
is the density and v the acoustical velocity of the excited mode in the 
medium, the impedance of Z at resonance becomes 


Z(nm) = 3(Z, + 23), (8) 
1 = koi and Zo = Rpsvr. 


The quantity & is given by 

4nd\? 1 
— aK) 3 ie) 
Since pv is the acoustical impedance of a medium, & is the proportion- 
ality factor between electrical and acoustical impedance. The quan- 
tity k depends entirely on the characteristics of the crystal alone, i.e. 
its dimensions and its electric and piezoelectric properties. For an 


X-cut quartz crystal 


K = 4.58, 
h = 14.3 X 10‘ c.g.s. e.s.u. 


k 


and 


: 2.86 X 108 
din mm, = eae a : 


where f is the fundamental resonant frequency. For a 30 mc./sec. 
crystal with an area of 1 cm?. we have 


Z, = 0.030(p1v1) ohms for (pv in c.g.s.). 


For mercury pv is 19.8 X 105 in c.g.s. so that each quartz mercury inter- 
face contributes 15,000 ohms to the impedance of the transducer at 
resonance. It is apparent from the form of equation (8) that at reso- 
nance the two bounding media act as series elements. 

From the direct approach employed in the preceding paragraphs, an 
expression for the impedance of the transducer, looking in from the 
electrical side, has been obtained. Its behavior at resonance has been 
investigated. A proportionality factor between electrical and acoustical 
impedance has been evolved and some actual numbers inserted. Also, 
expressions have been derived for the displacement amplitudes at the 
interfaces (equation 6a and 6b) in terms of the free charge on the 
crystal. A somewhat more general approach to the fundamental 
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equations will now be employed with a view to obtaining a compact 
equivalent circuit from which any specific result could be derived. 

One starts from equations (3) using parts (a) and (b) to eliminate 
Band C from parts c and d. The new equations are 


wpitiA = jwpv(A cot Bd — A’e~ #24 csc Bd) — - I, (9a) 


wp22A’e~ 1624 = jwov(—A csc Bd + A’e~ #24 cot Bd) + “s I. (9b) 


Next one identifies the displacements and pressures at the interfaces 
with fictitious voltages and currents respectively, as follows: 


A = —e,/h, A'e-ivxt = e,/h, 


e 5 ss h., 
Wp 1V,A = os) WPW2A e iBed = ait to. (10) 


The substitution of these quantities into 9a, 9b and 5 gives the following 
three equations: 
h? . 
jets (4; — I) = e, cot (Bd) + e2 csc (Bd), 
Piesg} ax 7 d 1 
ja*poS (2 — I) = e: csc (Bd) + e2 cot (8d), (11) 
4ndI 


V=ertes T hos 

These equations give the relations between the voltages and currents in 
a six-terminal network which could serve as an exact equivalent circuit 
for the transducer. As before, it will be expedient to change from a 
circuit where the electrostatic condenser is in series to one in which the 
condenser is in parallel and to neglect terms of the order of (wcZ)~, 
which is permissible when the crystal is damped. In the equations 
for the new network the roles of current and voltage in equation (11) 
are interchanged at each terminal and the impedance factors, such as 
h?/w*pvS, are multiplied by (&S/4md)? and become admittances. Con- 
sequently, the new network equations are 


(e: — V) = jZo[41 cot (Bd) + 2: csc (6d) J, 
(e2 — V) = jZo[t1 csc (Bd) + 22 cot (8d) |, (12) 
jukS 


ati Etats 2 ’ 


where Zy = kpv, the equivalent electrical impedance of the quartz. 
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Equations (12) can be represented by the equivalent circuit shown 
in Fig. 4, which first appeared in a report by H. Grayson.? Here the 
crystal is represented as a section of a transmission line (see the two 
parallel heavy lines) of electrical length 8d and characteristic impedance 
Z». By tying on impedances Z, and Z:, to the appropriate terminals, 
one can now rederive the specific results obtained in the first part of this 
section. 

The next problem is the investigation of the band-pass character- 
istics of the crystal. From the analytical expression for Z (6 and 7) one 
obtains the following formula for conductance, 


= 1\ _ {20 (cos (6d) — 1)* + ZZ; sin? Bd}(Z1 + Zs) | 
G = Re (3) oe Z(Zi1+ Z:2)? cos? (Bd) + (2:22 7 Ze) ain Bd |? (13) 


This quantity (13) is proportional to the power (V°G) absorbed by one 
crystal from a constant voltage source. The frequency response is 
obtained by varying fd. Figure 5 shows the behavior of (13) under 
four different conditions. All curves are periodic with period 27 for fd. 
When the crystal is acoustically matched on both sides, the curve is 
sinusoidal (A). The curve is narrowed somewhat (B) as the result of 
allowing one surface to vibrate freely (Z, = 0). If one employs a 


material of higher acoustical impedance on both sides (Z; = Z: = v2Z,), | 


a broader curve is obtained (C). Curve D is obtained when one of these 
surfaces is allowed to vibrate freely (Z. = 0) and shows a zero curvature 
at resonance, nearly the case of quartz with mercury on one side. This 
final result is somewhat unexpected and is of practical interest. 

.The expression (13) gives the total power absorbed by the crystal. 
For the two special cases considered above (Z; = Z2, and Zz = 0) 
there is no question how the power is divided between the two media. 
Though these are two situations of main practical concern, it is of 
interest to determine how the power is divided in the general case. 
To do this, we must return to equation (6) to evaluate the ratio 


Ofo| | déa| 
° ax |" Vo) & ax | 


vile 


This turns out to be 


Power transmitted to medium 1 
Power transmitted to medium 2 
By Z:2? sin? (Bd) oe ZZ [cos (Bd) net LP. (14) 
ZZ: sin? (Bd) + ZZ: [cos (Bd) was 1} 


It follows that the power absorbed by the medium 1 from a voltage 


? “The Quartz Crystal as a Generator and Receiver of Supersonic Vibration,’”’ H. Grayson, 
TRE Report T-1358. 
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Fic. 4. Transmission-line equivalent circuit. 


Pv. = PV pp Vo 

Pc = PV P, °° 
pve ya P= Pp, Vv, 
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Fic. 5. Acoustic output power of a crystal as a function of frequency. The crystal is 
assumed to be subject to a constant voltage. Fig. 1 gives meaning of subscripts. Resonant 
frequency of free crystal occurs for Bd = x. ‘“Q” refers to curvature at resonance. 
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source is proportional to 


ZZ; [cos (Bd) — 1} + Z,Z;? sin? (6d) 
[Z0(Z1 + Zs) cos (8d) }* + [(2:2Z2 + Zo) sin (6d) }? 
There is a similar expression for G: so that G; + G; = G. 


As an example the case of a quartz crystal with mercury on one side 
and water on the other has been worked out. In Fig. 7 the power 


G, (15) 


Ww 
” 
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oO 
a 
Mt) 
WW 
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at 


WwW tw 47 
6 6 63 


FREQUENCY (49) 
Fic. 6. Output power of crystal when backed by high impedance material. 


absorbed by the mercury, by the water, and the sum of these two quan- 
tities are plotted on a relative logarithmic scale as a function of 
frequency. 

In some applications the behavior of the frequency response in the 
neighborhood of resonance is of particular interest. One cannot use a 
simple series resonant circuit to represent both the amplitude and 
phase characteristics accurately here even over a limited frequency 
range. However, it is convenient to introduce a parameter, ‘“Q,” which 
is the Q value for the simple series resonant circuit whose amplitude 
response matches that of the transducer in curvature at resonance. 


II 
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For all but the lowest values, the parameter ‘‘Q”’ corresponds satis- 
factorily with the usual ratio of mid-frequency to distance between 
half power points. The formula for ‘‘Q”’ is given below and applies to 


(15) 


ne side 
power 


> 
3 
Pry 
~ 
e 
E 
° 
7 
» 
© 
on 
= 
E 
* 
c 
° 
ue 
~ 
A 
© 
é 


Power transmitted to water 


power 


Relative 


yuan- 
yn of , . ; ‘ 
jr r Sw iw lw 
a the Fic. 7. Frequency response; quartz crystal loaded with 
use a mercury on one side, water on the other. 
and 
ency the curve for power absorbed by the transmitting medium when the 
rhich transmitter crystal is driven from a low impedaance. 
ag ugn = meee. + Ze)? — (Z1 + Z:)(Zs* + 220) (16) 
a Z0(Z1 + Z:)? , 


ance. 


12 HuNTINGTON, EMSLIE AND HUGHES. {J. F. 1. 


where medium 1 is the transmitting medium. Here x is the order of the 
harmonic frequency of the crystal. For spetial cases the following 


simplifications occur: 
For Zi => + 


(17a) 
and for Z, = 0, 


“Q” = a Pez eZ.) }s (17b) 
1 


From (17b) it is evident for one free surface the response becomes 
antiresonant at frequencies Bd = nz (n odd), if Z1; > v2Zo. 

An interesting extreme occurs when one adjacent medium has a 
very large acoustical impedance (Z,-—> ©). The resonant points then 
occur at 8d = nz/2 for n odd, with a response characteristic (Fig. 6). 


a Zi sin’ (Bd) 
wee Z,? cos? (8d) + Z;? sin? (Bd) (18) 


It is apparent that the acoustical impedance of the piezoelectric 
material affects the bandwidth even though it does not enter into the 
expression for the equivalent electrical impedance at resonance (8). 
For quartz in mercury the response is broad. However, consider a 
10 mc. crystal driven at 30 mc. in contact with water on each side. 


pv for water = 1.5 X 10° c.g.s., 
pv for quartz = 15 X 10° c.g.s., 


“0” = a X 10= 23.5. 


Under the assumption of low-impedance electrical circuits used in 
conjunction with the crystals the shunt capacity of the crystals does 
not enter directly into the analysis of the preceding paragraphs. This 
assumption is not as drastic as might at first appear. Actually for 
many applications, particularly for delay lines where a broad and flat 
passband is desired, the crystal capacity is an element in a low-Q shunt 
resonant circuit. The impedance of this circuit is kept low in compari- 
son with Z,, Z., etc. by the bandwidth requirements. Under these 
circumstances the results of the last paragraphs are quite valid. 

It remains to calculate the actual loss through mismatch that occurs 
where a signal is transmitted through a delay line as depicted in Fig. 8. 
The equivalent circuit is given in Fig. 9. The load resistors shunting 
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the crystals at the transmitting and receiving ends are R, and R,, re- 
spectively. Consider the crystals’ Capacities as tuned. out. (The 
resolution of Z into two series elements }Z,; and }Z,,'is"possible only at 
resonance. ) 


If a current pulse of magnitude J is applied to the transmitting 
crystal, the voltage which is propagated down the line is given’by 
Zi 
Z1 + 22+ 4R. 


Here the transit time of the line is long compared with the pulse dura- 
tion, so that the pulse sees the characteristic line impedance. 


TRANSMITTER RECEIVER 


V = ZR, (19) 


ps PY WP’ 


Fic. 8. Delay line schematic. 


3, 
4 


Fic. 9. Equivalent circuit of delay line at resonance. 


If the other end of the line were terminated in its characteristic 
impedance, Z,/4 then V would appear across it. Therefore, looking 
back at the line, one sees a pulse generator of 2 V in series with the line 
impedance.* The voltage V, across the resistor R, is given by 


R, 
(Zi + Z2) + R, 


= IR; 


V, =2V 


8Z,R, 
(Z1 + Z: + 4R,)(Zi + Z: + 4R,) 


3 If the crystals at the two ends of the line differ with respect to any of the quantities upon 
which k (eq. 8a) depends, then Z; will have different values at the two ends, say Z,' and Z,’’. 
If we require that the available power be transformed without reflection by the line, then equa- 
tion (20) is valid (with appropriate superscripts in the denominator) if Z, in the numerator is 
replaced by the geometric mean of Z,' and Z,”. 


(20) 
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As has been remarked before, R; and R, can generally be neglected in 
comparison with Z, and Z,. The ratio of output to input voltage can 
then be simplified accordingly 

8Z:R, _ 
(Z: + Z:)? 


The corresponding power ratio is to the same degree of approximation, 


R.R,(8Z1)? 
(21 + Z:)' 


For R, = R, the power loss equals the voltage loss. 

To avoid multiple reflections in short lines, it is customary to try to 
make Z, equal Z; and to employ end assemblies which absorb the inci- 
dent energy. Equation (21) indicates that each such end assembly 
introduces a loss of 6 db. as compared to the situation where Z, is set 
equal to zero and one obtains 


(21) 


V,/IR, = 


; eer 2 | 

Where multiple reflections are absorbed by a sufficiently long path in 
mercury, or can be tolerated for other reasons, the loss in the line can be 
reduced by providing a low acoustic impedance, such as air, behind the 
crystals. The present practice is to back up the crystals with a solid 
electrode for mechanical support and trust that the contact is sufficiently 
rough to allow a layer of air at the interface. The predicted gain of 
6 db. per end assembly in signal strength in going from a completely 
absorbing to a completely reflecting end assembly has been experi- 
mentally observed. Also, for intermediate situations where partial 
reflection at an imperfectly absorbing back-surface, experimental re- 
sults were consistent with equation (21). 

The mismatch loss predicted by equation (22) has been experi- 
mentally verified for a tube of nearly one inch diameter. For smaller 
tubes the loss was greater than would be expected and might be due to 
mechanical edge effects at the periphery of the crystals. With a tube 
of 3’ ID the measured loss was about 10 db. more than predicted. 
With a tube of 25/64” ID the measured loss was found to be 6 db. more 
than predicted. 

The formulae developed here can be used to help choose optimum 
values for certain parameters. Two applications will be considered 
here, the choice of active crystal area and choice of fundamental or of 
harmonic frequency for the crystal operation. One desires to obtain 
the maximum ratio of output voltage to input current for constant 
bandwidth. In both cases the capacity at each crystal will be composed 
of two parts, a stray capacity C, which should be kept to a minimum and 
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a live capacity C,; = kS/4mnd, where S is the area of the crystal actively 

excited. The values of R, and R, must be inversely proportional to 

(C, + C,) if the bandwidth is to be maintained constant. The quan- 
tity we wish to maximize has already been shown to be 

V, 2s 8Z.R,R; 

T~G@tZ, (21a) 

For the determination of optimum crystal area S we note that Z; 

and Z, vary inversely as S (Eq. 8a). Under the conditions of varying 

S and holding d constant, the Z’s vary inversely as C;. It follows that 


a 

I (Ce +C,)?’ 
which as a maximum for C; = C,. This means that one chooses the 
crystal area so that live capacity matches the stray, for optimum ratio 
of voltage at receiver to current at transmitter. For optimum signal- 
to-noise the line capacity should be somewhat larger, depending on the 
characteristics of the amplifier’s first stage. 

Since, for constant S, the C; decreases with the resonant frequency 
of the crystal, excessive stray capacity is more costly at lower frequen- 
cies: The reason for this paradox is that the lower the frequency the 
more difficult it is to make the line capacity match a fixed stray capacity. 

To decide whether fundamental or harmonic operation is preferable 
one varies the crystal thickness d holding S constant. The Z’s vary 
directly as d? (Eq. 8a) or inversely as C;*. It follows that under these 
conditions 


(21b) 


V, Cp 
Lo 6. ee (ale) 


This means one wants the largest live capacity possible in this case, in 
other words, operation of the crystal on the fundamental frequency. 
For a given frequency the crystal that operates on the fundamental is 
the thinnest crystal that can be used. It is of interest to note that, it the 
stray capacity were negligible, working on the harmonics would corre- 
spondingly involve no additional loss. 

In the preceding analysis the driving and receiving circuits have been 
represented simply by a tuning coil and loading resistor in parallel. In 
any actual application the loss through mismatch can be reduced by a 
transforming network. Because of the bandwidth requirements, how- 
ever, it is often impractical to use as high a transformation ratio as 
would be required for match and the assumption that Z;, Z:2, etc. are 
large compared with the transformed impedance can be made. Up to 
this time it has not appeared feasible to use an acoustical transformer, 
particularly in the range of several megacycles. 
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C. WAVE PROPAGATION. 


It is of considerable importance to determine in as great detail as 
possible the way in which energy travels after leaving the transmitter 
of a supersonic system. To this end some calculations are first in- 
cluded on the radiation pattern to be expected from a piezoelectric 
crystal transmitting single-frequency compressional vibrations into a 
semi-infinite medium. Secondly, we shall make some remarks on the 
wave pattern present in a semi-infinite cylindrical tube at the end of 
which a similar crystal is vibrating. The second case corresponds more 
closely to the actual physical set-up for a delay line, though results ob- 
tained from the first case are also of value in determining mechanical 
tolerances for delay line construction. 

Two assumptions in particular which are made in connection with 
both cases should be noted. First, as is customary in such calculations, 


P, (Xo,Yo, Ze,) 


Fic. 10. Geometry of the radiation problem. 


we determine the beam pattern in a certain region of “‘empty’”’ space 
and then assume that the wave function will remain unchanged if the 
receiver is placed in this region. This is of course valid for pulse trans- 
mission which was our main concern at the Radiation Laboratory. 
Secondly, we consider that out medium is non-viscous. In the next 
section of this part, certain effects of viscosity will be considered. 


Case I: The Radiation Pattern in a Semi-infinite Medium. 
General and Fresnel Case 


The nomenclature for the problem is indicated in Fig. 10. A 
crystal with a circular cross section of radius p; and with its center at 
(0, 0, 0) is executing small amplitude, single-frequency vibrations in the 
direction of the x-axis. The portion of the YZ plane for which p> p: is 
assumed to be a rigid baffle. We wish to know the radiation pattern as 
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observed at an arbitrary point (Xo, Yo, Zo) to the right of the YZ plane. 
It is well known that the velocity potential, W, for sound waves of 
small amplitude in a perfect compressible fluid satisfies the scalar wave 
equation.‘ 
ew 
of 


where v is the velocity. The boundary conditions on the plane X = 0 
are: 


= vv?W, (22a) 


as constant # O(p < po) 
—> = 22b 
> ae Eee (o> pn): - 


These boundary conditions neglect the fact that the face of the crystal 
does not always lie in the same plane. But we shall always calculate 
the amplitude of the sound wave at a distance from the circular opening 
which is large compared with the wavelength, and in this case this is a 
valid approximation. We use the Kirchhoff condition which states: ® 


ow 
t 


‘on y 
Vo = fe 2 — sin 2r 77 x) cos (ur)deo, (23) 


where, referring to Fig. 10, 


Vo = velocity at Po, the observation point, 
ow ?A., 3 
—— ='— in which A is a constant, 
On r 
\ ="wavelength, 
T =|period. 
(nr) =fangle between direction of X-axis and r, 
t = time variable, 
do ="element of area of surface x = 0. 


Referring to Fig. 10 we see 


2 2 a _ 
ponte ares 2ZZ0 


By analogy to light diffraction problems we distinguish two cases. 
If it is assumed that 7 is expanded as 


2 YV¥o+Z2 
r=R(14+55 e = ‘). 


4A. G. Webster, “Dynamics”; Stechert, pp. 542-3 (1922). 

5 For discussion and partial derivation see Stratton’s “Electromagnetic Theory,’’ McGraw- 
Hill, pp. 460-463 (1941). It should be pointed out that our radiation problem resembles that 
of determining the diffraction pattern to the right of a circular opening upon which a plane 
wave of light is incident from the left. 
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we are dealing with the Fresnel case. If 
r=R(1 es Eye 22s) 


we are dealing with the Fraunhofer case. We consider first the Fresnel! 
case. Hence (23) can be simplified to (24). 


ve=a' f sin2e (4-*) de, (24) 
opening : 


A cos (nR) | 
MEAS 


This integral is now in a form which has been evaluated. With the 
omission of the factor A’ the integral can be written as 


[Of sin ax 5-2) piddp = Msin(o+2rt), ( 


where p and ¢ are polar coordinates for the circular opening, M and ¢ 
are functions of Xo, Yo, and Zo. The elaborate and laborious process 
of computing M for several values of (Xo, Yo, Zo) has been carried 
through in a paper by E. Lommel on “The Diffraction Phenomena of a 
Circular Opening and a Circular Screen.’”’* These calculations check 
the experimental data on diffraction of light towithin about 1 or 2 percent. 


A’ = 


Fraunhofer Case 
The solution of the Fraunhofer problem is well known,’ 


Ye * V1 

1 ( Fo’) 
2s :) 
r pix 


where J; is a Bessel function of the first kind and x is the sine of the 
angle which R makes with the X axis. 


On the Axis 
An analytical expression for the disturbance on the X-axis for the 
Fresnel case can be obtained quite readily. By (23) for a plane wave 


incident normally. 
sin 24 (7 x) de do 
yam am (27) 


A opening 


(26) 


‘E. Lommel, Abhandlungen der Bayerischen Akademie der Wissenchaften, 15, 233 (1886). 
7M. Born, “Optik,” pp. 159-161 (J. Springer, 1933); or P. Morse, “Vibration and Sound,” 
p. 257 (McGraw-Hill, 1936). 
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This expression is valid so long asr > p;. It can be integrated directly 
with respect to 7 and gives 
d 


~ 2h eee he Sige 
Vs = 24 sin (22) sin 2x rx) (28a) 
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(26) 


f the 


* 2 3 ‘ 
X—» (Distance from plane of sound source in PAY 


Fic. 11. Intensity of sound on X-axis vs. distance from sound source. 


where A is as before the amplitude of the velocity at the crystal. This 
equation is plotted in Fig. 11. At large distances the amplitude of the 
velocity reduces to the well-known formula: 


Vo = ae (a pi’). (28b) 
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Total Force on Receiving Crystals 


Since the wave incident on our receiver crystal is not a plane wave, 
pressures of different amplitude and phase will fall at different points 
of the receiver crystal. It appears that the factor of importance in 
determining the response of the receiver crystal is the total force on it. 
This quantity F can be determined by integrating over the area of the 
receiving crystal. 


r=Pf Msin (0 +2" 4) ao, (29a) 


Fl = P| f M cos sie | +| f M sin ode J, (29b) 


where P is the pressure amplitude at the transmitting crystal. 
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Fic. 12. Force on receiver crystal (p2 = 1) vs. distance between crystals. The two 
crystals have equal radii and are aligned parallel. Dotted —-—-—- gives asymptotic be- 
haviour at large distances. 


Two cases are shown plotted in Figs. 12 and 13. At large distances 
the force on the receiver crystal approaches asymptotically to the 
hyperbola 


(30) 


which is also drawn in Figs. 12 and 13 (ps2 is radius of receiver crystal). 
Fig. 12 applies to the usual case of two crystals of the same size, or one 
crystal and a plane reflector. For Fig. 13 the radius of the receiver 


crystal is 4 that of the transmitter. The curve shows an oscilla- 
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tory motion at small separation and a pronounced maximum near 
d = 1.2p,2/X, at which point presumably it is being effectively irradiated 
by all parts of the transmitter crystal. The value of the average pres- 
sure here is nearly 60 per cent. greater than the value, P, at the trans- 
mitter. 


Case II: The Radiation Pattern in a Semi-Infinite Tube with a Vibrating Crystal at One End. 


The analysis of ultrasonic propagation in a tube of semi-infinite 
length depends on the choice of boundary conditions at the wall of the 
tube which correspond to physical reality. If the boundary condition 
taken is that the normal component of the velocity at the wall of the 
tube shall be zero, then for a system with a transmitter equal to the 
tube radius a nearly plane wave can be made to travel down the tube. 
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Fic. 13, Force on receiver crystal (p2 = 491) vs. distance between crystals. 
The radius of the receiving crystal is one-third of the transmitter. 


If the radius of the transmitter is less than the radius of the tube, the 
disturbance in the tube can be considered as the sum of the disturbances 
due to an infinite number of cylindrical modes travelling down the tube 
with differing velocities. 

Another boundary condition which may be reasonable in some cases 
is that the pressure at the wall of the tube is zero. Such a boundary 
condition may be approximately valid when there are air spaces between 
the fluid and the tube (as often occur for mercury in a steel tube). 
With this boundary condition it is found that a series of cylindrical 
modes will travel down the tube with varying velocities. We are 
indebted to Mr. H. J. McSkimin of the Bell Telephone Laboratories 
for allowing us to see his unpublished analysis of this situation. His 
results show that the intensity of the ultrasonic beam falls off very 
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slightly at large distances in the tube in contrast to the inverse square 
law spreading of the free-space beam. 

In Section D where viscosity effects also are considered there will be 
presented some experimental data on the question of propagation of 
waves in tubes. 


D. ATTENUATION. 
Free Space Beam 


A “free space” beam in a liquid decreases in two distinct ways. 
There is first the spreading of the beam due to diffraction, which gives 
rise to an inverse square law for the intensity as large distances from the 
transmitting crystal (see eq. 28b of preceding section). The amplitude 
behavior can be described by 

a=a-7, (31) 
where 4» is the amplitude near the transmitting crystal, S the area of the 
crystal, d the distance from the crystal and \ the wavelength. Eq. (31) 


holds provided dx > S. 
The second kind of decrease is an absorption due to viscosity and 


thermal conduction in the liquid, which results in an exponential law 


for amplitude versus distance of the form a = a,ye~**. The theoretical 
formula * for a is given by the expression. 


io asf * 2%. SUA |. (32) 


pv? L3 c. 
where 7 is the viscosity coefficient, f the frequency, p the density and v 
the acoustic velocity. In the second term y is the ratio of the specific 
heats, K the thermal conductivity and c, the specific heat at constant 
pressure. For non-metallic liquids the effect of thermal conduction 
is negligible. For many liquids this simple theory predicts an attenua- 
tion much smaller than that found by experiment. For mercury, 
however, which was widely used in delay line operation, agreement is 
good. Theory predicts an attenuation of about 1.1 db. per ft. at 10 mc., 
whereas the experimental values indicate 1.2 + .2 db. per ft. Here 
the heat conduction loss accounts for 80 per cent. of the attenuation. 
Since the viscosity coefficient decreases with rising temperature, 
attenuation generally decreases also. The effect is very marked in the 
case of water.°® 


8 W. P. Mason, “‘Electromechanical Transducers and Wave Filters,” Van Nostrand, p. 305 


(1942). 
°C. E. Teeter, Jr., J. Acous. Soc. Am., 18, 488 (1946). 
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Guided Beam 

It seems to be possible to avoid at least in part the decrease due to 
spreading of the beam by confining the beam within a tube of the same 
diameter as the transmitting crystal. This is what is meant by a guided 
beam. If the liquid could slip freely over the tube surface, a piston- 
like motion would take place corresponding to perfect plane-wave 
propagation, without any loss due to spreading. Actually, the condi- 
tions at the wall are not as simple as this and depend on the nature of 
the surface. A possible assumption is that the liquid sticks to the sur- 
face instead of sliding freely over it. Crandall }!° has shown that the 
viscous losses take place in a very thin layer, of the order of 10~-* cm. 
in thickness, across which a major portion of the velocity change occurs. 
The thickness of the layer is given by 


$= (.)' (33) 


where 7 is the viscosity, p the density, f the frequency. The drag 
produced by the layer on the liquid column causes an exponential at- 
tenuation with an amplitude absorption coefficient per unit length given 


by 


p= os (34) 


br ’ 


where 6 is the radius of tube and \ the wavelength. For 10 mc. waves 
in mercury contained in a tube of radius 0.3 cm. we have 6 = 2.85 X 107° 
cm. and 6 = 0.0040 cm. Some experiments have been carried out 
at 10 mc. which indicate that the attenuation of ultrasonic wawes in 
mercury shows a dependence on tube diameter of this order of magni- 
tude. The effect was, however, influenced by the character of the 
tube’s inner surface. It has already been mentioned in Section C 
that the boundary conditions may have a profound effect. For ex- 
ample, if, instead of a high acoustic impedance in the steel wall, the air 
occluded in the irregularities of wall surface offers a low impedance, 
there will be no energy absorption at the tube wall and the propagation 
will always consist of many cylindrical waves travelling with different 
velocities down the tube. 

Since the guided beam is not subject to the attenuation from spread- 
ing given by equation (31), total attenuation is generally less in this 
case than for the free space beam. Only for long distances can the 
tubular attenuation give an effect comparable with the spreading. Of 
course both beams are subject to the absorption varying as the square 
of the frequency (eq. 32). 


107. B. Crandall, ‘““Theory of Vibrating Systems and Sound,” Van Nostrand, p. 229 (1927). 
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Science Captures a Chinese Art.—In China the making of soya sauce is 
mainly a household art, with grandmothers passing along to daughters and 
granddaughters the details of recipes and the strains of ferments used in the 
process. 

In the United States, soya sauce—as is true of many other food products— 
is chiefly a factory product. Each manufacturer is concerned not only with 
producing a sauce of high quality, but also a standardized sauce in which each 
batch will have the flavor, color, and taste of every other batch. This calls 
for standardization of the raw materials and for use of standard strains of 
microorganisms in the fermentation of the sauce. 

The U. S. Department of Agriculture has announced that strains of four 
organisms desirable in preparing soya sauce have been added to the culture 
collection of industrial ferments at the Northern Regional Research Laboratory 
at Peoria, Ill. These include two molds, a yeast, and a bacterium. They will 
be maintained as pure cultures and will be available to industrial users. This 
makes it possible for a fermenter to make a fresh start with new and pure cul- 
tures, if at any time his stock cultures become contaminated with ‘‘wild’’ molds 
or yeasts that injure the quality or uniformity of his product. 

The Fermentation Division of the Bureau of Agricultural and Industrial 
Chemistry credits Mr. Pei Sung King of the National Bureau of Industrial 
Research, Chungking, China, with aid in selecting the strains of organisms and 
standardizing a process of fermentation that yields a high-quality sauce. Mr. 
King has been a guest worker at the Northern Laboratory, and suggested 
various methods of fermentation which the scientists tested and compared in 
working out the method they now recommend to manufacturers. It is not 
practical for home use or for small quantities of sauce. 

The preparation of soya sauce calls for a brine fermentation of the beans for 
from 30 to 90 days. But the “starter’’ used in this process is a mixture of three 
previously prepared cultures; (1) of a mold that develops on cooked rice; (2) of 
a yeast working on soya broth; and (3) of a bacterial fermentation of soya 
broth. For a quality product these should be mixed at the right stage of 


development in suitable proportions. 
Ri HH. O. 
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AN INTEGRAL-EQUATION APPROACH TO PROBLEMS 
OF VIBRATING BEAMS.* 


BY 
WALTER T. WHITE, Sc.D.} 


PART I. 
ABSTRACT. 


In this paper integral equations are applied for the calculation of the normal modes of 
vibrating beams. Both exact and approximate methods of solving the integral equation are 
considered. The Green’s function, or kernel, of the integral equation is constructed for both 
uniform and nonuniform beams. Solutions for the normal modes of a uniform cantilever are 
given. A nonuniform, naturally-twisted turbine blade is studied in detail and the first and 
second normal modes are calculated by the integral-equation method. 


1. INTRODUCTION. 


In many mechanical designs, it is important to eliminate vibrations 
associated with mechanical resonance. Calculations of the normal 
modes of the system often reveal objectionable resonances, and when 
they are known, methods can be devised to reduce or eliminate them. 
In some cases, damping or appropriate shock mounting may achieve the 


desired result. But in others, it is necessary to change the inertia and 
stiffness of the system. The theory of vibrating beams, when appli- 
cable, is the most effective basis for such calculations. 

Vibrating members treated as beams are rotating shafts, turbine 
blades, hulls of ships, airplanes, missiles, and other complex structures. 
The beam and its supports may be sufficiently simple for straightforward 
analysis by classical methods with differential equations. If the prob- 
lem is complicated by variations of physical properties, the solution by 
integral equations gives results more conveniently than solution by 
differential equations. The formulation of an integral equation is 
based upon a simple picture of the beam and the application of the 
superposition principle. Variations of physical properties and such 
effects as rotary inertia, longitudinal inertia, and gyroscopic moment 
cause no particular difficulty with the integral-equation formulation. 

The purpose of this paper is to outline the integral-equation ap- 
proach to the calculation of normal modes of vibrating beams. Certain 
properties of integral equations required for the solution of practical 
problems are given. A method is developed for constructing the kernels 
for the integral equations of vibrating beams. To demonstrate the 
integral-equation method and its accuracy, an application is made to a 

* This paper is from a thesis, “Integral-Equation and Approximate Solutions for Normal 
Modes of Vibration,”’ submitted to the department of Electrical Engineering, M.I.T., 1941. 

t Research Engineer, Sperry Gyroscope Co., Great Neck, N. Y. 
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simple cantilever. A second application is made to a nonuniform 
naturally-twisted turbine blade. 


2. FORMULATION OF THE INTEGRAL EQUATION. 


Two procedures generally are available for the integral-equation 
formulation of a vibration problem: (1) obtain the differential equation, 
and then make a transformation to the integral equation; (2) formu- 
late the integral equation directly from the physical nature of the 
problem. Evidently the second procedure has distinct advantages. 
Indeed, one of the principal advantages of the integral-equation treat- 
ment is the closeness with which it is tied in with the physical picture. 


y 
j Br 9 wt 
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+f 


Fic. 1. Inertial loading on vibrating cantilever. 


As an example, consider the cantilever. To formulate the integral 
equation, it is necessary to consider the inertial loading while the canti- 
lever is freely vibrating. If the cantilever is observed at the instant of 
its maximum downward deflection, as in Fig. 1a, each elemental mass 
has a deflection ++(z,#) and an upward acceleration —a(z,#). The 
product of mass m per unit length and acceleration gives the upward 
inertial force which the beam exerts upon its mass. By the principle 
of action and reaction, the mass exerts a downward force or loading 


f(z, t) per unit length that is given by 
f(z, t) = — ma(z, 2). (1) 
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The loading on an elemental length Av is equivalent to a concen- 
trated load f(v)Av (see Fig. 1b). From texts on strength of materials,*:! 
the static deflection of a uniform cantilever caused by a unit-concen- 


trated load applied at z = v is 


G(z, v) = (30 — 2) /6EI, 20 (2) 
= v(32 —v)/6EI, v<32}’ 


where EJ is the flexural rigidity of the cantilever. The static deflec- 
tion caused by a unit-concentrated load generally is called a Green’s 
function. The deflection caused by the loading fAv is 

Ay(z, v) = G(z, v)f(v) Av. (3) 


By superposition, the deflection at z is given as the summation of the 
deflections caused by each elemental loading along the beam, or 


we. 1) = f " Gls, vd fle, Dae. (4) 


The time variable ¢ is introduced, since the loading f in Eq. 1 is a time- 


varying quantity. 
Since undamped oscillations of small amplitude are harmonic func- 
tions of time, y(z, ¢) can be written as 


y(z, t) = y(z) sin wt, (5) 


where w is the frequency of oscillation in radians per second. The ac- 
celeration is obtained by differentiation of Eq. 5, and is 


a(z, t) = — w*y(z) sin of. (6) 
Substitutions from Eqs. 1, 5 and 6 in Eq. 4 gives 


y(z) sin wt = mw* r G(z, v)y(v) sin widv. (7) 


Since the time function is not dependent upon the variable of integra- 
tion, it is cancelled to give 


y(z) =X f £ G(z, v)y(v)dz, (7a) 


where the constants are lumped in the parameter \. 

Equation 7a is the integral equation for the freely vibrating canti- 
lever. The kernel, or Green’s function, G(z, v), is given by Eq. 2, and 
y and \ are unknowns. The problem is to find values of \ for which 
non-zero functions y exist. Equation 7a is called a homogeneous integral 
equation; the d’s are called eigenvalues and the corresponding y’s are 
called eigenfunctions. The eigenvalues and eigenfunctions determine 
the normal modes of vibration for the cantilever. 


* For numbered references see Bibliography. 


Water T. WHITE. 


3. PROPERTIES OF HOMOGENEOUS INTEGRAL EQUATIONS. 


Before completing the solution of Eq. 7a for the cantilever, it is well 
to consider certain properties of homogeneous integral equations’. 
which are useful in the solution of practical problems. 

A general form of the homogeneous integral equation is 


y(s) =X 7 K(z, v)y(v)dv, (8) 


where a and 6 are constants and K(z, v) is the kernel corresponding to 
the Green’s function in Eq. 7a. In vibrations and related problems, the 
kernel is often a symmetric function of z and v: 


K(z,v) = K(v, 2). (9) 


If the kernel is asymmetric, it can be made symmetric by a simple 
transformation, as will be shown later. In practical problems the kernel 
is real, continuous and not zero on the interval (a, 0). 

An obvious solution of Eq. 8 is 


y(z) = 


Although this result yields no useful information, useful solutions can 
be found. For certain eigenvalues \ there exist corresponding solu- 
tions or eigenfunctions y which are written 


yi(s) = ry i: K(sz, v)yi(v)dv 


As af at Vv) yo(v)dv 


=), ff K(z, v)y,(v)dv 


| 
} 
) 


Clearly, each eigenfunction is determined only within a constant multi- 
plying factor, for, if y,(z) is a solution of Eq. 8, substitution shows that 
cy,(z) is a solution where ¢ is a constant. Solutions in Eqs. 10, there- 
fore, only determine the form of the eigenfunctions. The eigenvalues, 
however, are uniquely determined. 

A general property for eigenfunctions of symmetric kernels is their 
orthogonality. This property is expressed as 


[- ya(o)ya(o)d0 = 0, (m # n). (11) 


This s 
all vah 
bilinea 
in tert 
the un 
Oft 
which 
is asso 
of the 
cantile 


and th 


Obviou 


Byas 
metric. 


andar 


With t 


Since t 
integra 
tion fre 


1g to 
, the 
(9) 


mple 
-rne] 


Jan., 1948.] PROBLEMS OF VIBRATING BEAMS. 


Ifm = n, the integral has a value NV, called the norm; that is, 


f * Cya(o) Plo = N,. (12) 


A useful property of the kernel is its bilinear expansion. For an 
infinity of eigenvalues and eigenfunctions, the bilinear expansion is 


given as 


K(z, 0) = 5 y.(s)y.(0)/r.N.. (13) 


n=) 


This series expansion for the kernel must be uniformly convergent for 
all values of z and v on the interval (a, 0). The expansion is Hilbert’s 
bilinear formula.? It gives the kernel in the manner of a Fourier series 
in terms of the eigenfunctions. When applied to physical problems 
the uniform convergence is assured. 

Often the kernel is asymmetric. For example, if the cantilever 
which led to Eq. 7a has a nonuniform cross section and loading, there 
is associated with the Green’s function a function »(v) which takes care 
of the variations of mass and flexural rigidity along the length of the 
cantilever. The integral equation then has the form 


y(e) = 2 f° Gl, v)a(o)y(o)ar, (14) 


and the kernel is 
K(z,v) = G(z, v)u(v). (15) 


Obviously, the kernel is asymmetric, since 
K(z, v) ¥ K(v, 2). 


By a simple transformation, the asymmetric kernel can be made sym- 
metric. Multiply Eq. 14 by Vu(z) and define a new eigenfunction 


Y(z) = y(z)Vu(2), (16) 
and a new kernel 


K,(z, 0) = G(z, v)Vu(2)u(2). (17) 
With these substitutions Eq. 14 is written 
¥(i) «> f ” K.(a, v) ¥(v)dv. (18) 
0 
Since the kernel of Eq. 17 is symmetric, Eq. 18 can be treated as an 
integral equation with a symmetric kernel. The orthogonality condi- 


tion from Eqs. 11 and 16 is 


[o u@)ya(oya(0)do = 0, (m # 2) (19) 
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The foregoing properties are useful for the solution of practical 
integral equations. However, complex kernels arise in certain problems 
and such problems require special orthogondlity conditions and a more 
complicated bilinear expansion. After the formulation of the problem, 
an inspection will reveal whether the integral equation is similar to 
Eq. 8 or 14, or whether it is more complex and requires the derivation 
of a special orthogonality condition or bilinear expansion. 


4. SOLUTION OF INTEGRAL EQUATIONS. 


Practical methods for solving integral equations seldom are exact. 
More often solutions are obtained by approximations, and with reason- 
able effort such solutions can be calculated with an accuracy comparable 
to the physical data. Two methods frequently used for solving integral 
equations are: 


1) Obtain a solution by an intuitive guess, and prove it by sub- 
stitution. 
2) Obtain a solution by successive approximations. 


As an example of solution by the first method, consider the cantilever 
of Eq. 7a. Substitution of the Green’s function from Eq. 2 gives 


y(z) =X : v*(3z — v)y(v)dv + d 2 2?(3v — 2)y(v)dv, 


where 
= mw*/6EI. 


As a guess for the solution, try the equation 
y(z) = Boosh kz + Csinh kz + Dcos kz + E sin kz. 
Substitution in the integral equation gives the eigenvalue 
= k*/6, 
and the eigenfunction 
y(z) = cosh kz — cos kz + 6 (sinh kz — sin kz). 
The constant } is 


b = — (sinh kL — sin kL)/(cosh RL + cos kL) 
= — (cosh kL + coskL)/(sinh RL + sin RL) }’ 


where kL is a root of 
cosh kL coskL + 1 = 0. 
A few values of kL are 
kL = 1.8751, 
koL = 4.6941, 
R3L = 7.8548. 


From the foregoing the first three natural frequencies of the canti- 
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lever are obtained from the eigenvalues as 


w, = (3.516/L2)VEI/m, 
w, = (22.03/L*)VEI/m, 
ws = (61.70/L*)VEI/m radians per second. 


Generally, a guess as in the foregoing will not lead to the solution in 
a closed form. One reason for this is that the Green’s function may be 
known only as a family of curves. For problems not easily solved in 
closed form, the method of successive approximations is appropriate. 
When used with numerical or machine integration, the method of 
successive approximations is applicable to problems in which the 
kernel is known from tabulated data or plotted curves. 

In the method of successive approximations, the normal modes are 
calculated singly beginning with the first, or lowest mode, and con- 
tinuing in order to the higher modes. For the first mode, an approxi- 
mation to the eigenfunction is used as a generating function. It is 
necessary only that the generating function be not orthogonal to the 
first eigenfunction. Although it is improbable that a function arbi- 
trarily selected would be orthogonal to the unknown eigenfunction, it 
can be shown that a Green’s function is certain of meeting the non- 
orthogonality condition. Accordingly, the generating function go(z) is 


taken as gidak mm Olea), (20) 


where @ is an arbitrary value which gives a Green’s function not identi- 
cally zero. 

For the first eigenfunction, select for a the value where judgment 
indicates a concentrated load to be most nearly equivalent to the dis- 
tributed inertial load. For the uniform cantilever the value is approxi- 
mately a = 0.8L. (The best value can be shown to be a = 0.8021L.) 
When the generating function is selected by this method, solutions are 
obtained with a small number of successive approximations. There- 
fore, unless some other function obviously has advantages, a Green’s 
function is used as the generating function. 

From the generating function go(z), a set of equations is formed by 
successive iterations with the kernel. For the cantilever the itera- 
tions follow from Eq. 14 as 


g:(z) = f * Gis, WulOeclolde 


) 


> 


ges) = f° Cle, v)ulo)gu(o)de 


gai(e) = f° Cle, v)a(o)gn(o)ao 
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As n becomes large, these successive iterations give the solution of the 
integral equation 2 as 


yi(z) = g(2). (22) 


This is the eigenfunction of the first mode. The first eigenvalue jis 
calculated by substitution of Eq. 22 in Eq. 14 which gives 


ga(s) =a J “Gls. oyaloletende: 


Substitution from Eqs. 21 gives 
Zn(Z) = Ai¥nsilzZ). 


The eigenvalue }, is calculated from Eq. 23a for any value of z. How- 
ever, where approximate methods are used, a more accurate value of 
d1 is calculated as the average value 


a i gn(v)dv a 
: fi ennioao es 


Usually the first eigenvalue is calculated after one to three iterations. 
The first eigenfunction generally is calculated after two to four itera- 
tions. Since the convergence is very rapid, the calculation of the eigen- 
function is complete when the form of g,:(z) is similar to g,(z). 

Calculation of the second eigenvalue and eigenfunction by approxi- 
mate methods is possible only after calculation of the first eigenfunction. 
A generating function is selected which is orthogonal to the first eigen- 
function and not orthogonal to the second eigenfunction. The first 
condition is met by a procedure to be developed below. As for the 
second condition, it is improbable that an arbibrary function would be 
orthogonal to the unknown eigenfunction; yet, where certainty is 
desired, a Green’s function given by Eq. 20 is selected. A good value 
to choose for a is the estimated coordinate where the maximum deflec- 
tion of the beam occurs when it is vibrating in the second mode. 

To construct a generating function go(z) orthogonal to the first 
eigenfunction y,(z), form the equation 4 


go(z) = fo(z) — ay,(z), (25) 


where @ is a constant and f(z) is an arbitrary function, preferably a 
Green’s function G(z,a). Multiplication of Eq. 25 by u(z)y.(z) and 
integration over the length of the beam give 


[i or @aete)as = fiulsyst@) fleas — @ f* nledyale)ds. (26) 


The integral on the left is the required orthogonality condition and is 


Ai 
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zero if 


J ” u(2)y1(2) fol2)de 
J ” u(z)y2(2)de 


Accordingly, the generating function for the second eigenfunction can 
be constructed from Eq. 25 with a evaluated by Eq. 27. 

Successive iterations of the generating function for the second eigen- 
function give 


(27) 


5 


L 
gi(z) = J G(z, v)u(v)go(v)dv 


« 


£nii(Z) = f G(s, ) u(v)ga(v)do | 


If each integration is exact, the mth iteration approaches the second 
eigenfunction as m becomes large. However, exactness is possible only 
with formal integration, and if approximate integration is used, a 
special procedure is required to assure the convergence of Eqs. 28 to the 
second eigenfunction. The procedure simply is to make certain that 
each successive generating function is orthogonal to the first eigen- 
function. This requires a construction similar to Eq. 25 after each 
iteration to obtain an orthogonalized generating function, g’,(z), given 


by 
£'n(2) = £n(2) — anyi(2). (29) 
The (x + 1) iteration gives 


L ; 
Zasi(z) = I G(z, v)u(v)g’ .(v)dv. (30) 
As n becomes large, the second eigenfunction is the orthogonalized 
g'n41(2) and is given by 
Y2(Z) = n4i(Z) — On4191(2). (31) 
The second eigenvalue is given by a relation similar to Eq. 24 as 


L 
f g’ .(v)dv 
eA 0 
- , 
J o' axi(o)dv 


Another method for calculation of the second eigenfunction makes 
use of the bilinear expansion, Eq. 13. Since Eq. 13 is for a symmetric 
kernel, the integral equation must be in the form of Eq. 18: 


Xe 


(32) 


V(z) = f ” Ka, v) ¥(v)do. (18) 
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The bilinear formula gives 


Oe 
K,(z, v) = > Y,(z) Y,(v)/A.N a. (13a) 

n=l 
If the first eigenfunction and eigenvalue are known, they are subtracted 

from K,(z, v) to give the perturbated kernel K,(z, v), 
K,(z,v) = K,(z,v) — Yi(z) Vi(v)/A.Ni. (33) 
Since the remaining eigenfunction components in K,(z, v) are orthog- 
onal to Y;(z), no error results in the calculation of the second eigenfunc- 
‘tion if the pertrubated kernel, K,(z,v), is used instead.of K,(z, v). 
Hence, with a suitable generating function, which need not be orthog- 
onal to the first eigenfunction, the second eigenfunction is calculated by 
successive iterations to give 


gs) = f° KC, o)go(o)de 


£asi(s) = [' KG v)g,(v)dv 


There is no tendency to converge to the first eigenfunction and succes- 
sive generating functions are not required to be made orthogonal to the 
first eigenfunction. For large values of m the second eigenfunction is 
given by 

Y;(z) “se g1(2), (35) 


and the second eigenvalue is 


(36) 


The extension of the foregoing method to the calculation of higher 
modes is straightforward. A new perturbated kernel is formed by 
subtraction of previously calculated components from the bilinear 
formula, and the next higher mode is calculated by relations similar to 
Eqs. 34, 35 and 36. Alternatively, the higher modes are calculated as 
in Eqs. 28 and 29 provided the generating function for each successive 
iteration is made orthogonal to the previously calculated lower modes. 

An application of the foregoing methods is made for the calculation 
of the first two normal modes of the uniform cantilever. Although the 
complete solution was obtained for the integral equation, the method of 
successive approximations is used for purposes of illustration. A com- 
parison of numerical results from the approximate method with the 
exact solution illustrates the accuracy of the method. 
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Several mearis are available for the numerical evaluation of the 
integrals. A standard calculating machine can be used with Simpson’s 
rule. Machines, such as the cinema integraph® or the differential 
analyzer,* if available, can be used for evaluation of the integrals. In 
this example the differential analyzer is used. 

The integral equation to be solved is 


y(2) = f ” K(s v)y(v)do, 
| 


where K(z, v) is a symmetric kernel and is obtained from Eq. 2 as 
K(z, v) = 2(3v — 2), S|. 
= y*(3z — v), v<2| 
\ is the eigenvalue which includes the unknown natural frequency and 
the constants of the equation. In this example 
\ = mw’/6EI. 

If the integral equation is non-dimensionalized by the substitutions 

z=zL 

vy=v'L}’ 
it is written in non-dimensional form (with the primes omitted for 
convenience) as 


y(z) = i ' K(e, 0)y(0)do, 


and 
dh = mL? /6EI. 


Numerical values for the kernel of the integral equation are in Table I. 
From the estimated form of the freely vibrating cantilever, the generat- 


TABLE I, 

Symmetric Kernel for Uniform Cantilever. 
K(z, v) = #(30 — z), 2 < v, K(z, v) = K(v, z) 
(Multiply by 10-*) 

0.2 0.3 0.4 0.5 0.6 0.7 0.8 


2 
° 


432 

540 686 

648 833 1024 

756 980 1216 1458 

864 1127 1408 1701 2000 
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ing function is taken as 
Zo(z) = kz?, 


where & is any convenient constant. Two successive iterations of the 
generating function, as in Eqs. 21, are sufficient to evaluate the eigen- 
function. The results from the second iteration are in Table II, as are 
the exact values. The eigenvalue is calculated from Eq. 24, and the 
natural frequency, calculated from the eigenvalue, is given in Table I]. 
The errors between the exact and calculated values for both the eigen- 
function and natural frequency are negligible. The normalization con- 
stant is calculated from Eq. 12 and listed in Table II. 

For the evaluation of the second mode, the perturbated kernel is 
calculated from Eq. 33 as 


K,(z, v) = (30 — 2) — yilz)ai(v)/AiMi, <0. 


The values of i(Z), Aa and JN, tabulated in Table II for the second 


TABLE II, 
Normal Modes of Uniform Cantilever. 


| 
1 (s): First Mode yz (s): Second Mode 
| 


Third 
Iteration 


Second 


Iteration Exact 


Exact 


0.0000 0.0000 0.000 
0.0168 0.0168 —0.092 
0.0639 0.0639 —0.301 
0.1365 0.1366 
0.2299 0.2302 
0.3395 0.3391 
0.4611 0.4612 
0.5908 0.5905 
0.7255 0.7247 
0.8624 0.8635 
1.0000 0.9997 1,000 


meessssssss 
SOBUIHDAURPwWHHROS 


w: Natural frequency in radians per second 
3.5160 VEI/m 3.5176 VEI/m 22.03 VEI/m 22.04 VEI/m 
L L L L? 
N: Normalization constant 
2.001 2.001 


iteration are used for calculation of K,(z,v). The generating function 


for the second mode is 
go(2) = kK,(z, 0.6), 


where¥k"is a convenient constant multiplier. The results from the 
third iteration of the generating function are compared with the exact 
values in Table II. The errors for the approximate calculation of the 
second mode are almost as small as the errors for the first mode. 


(To be continued in next issue.) 
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VISUAL MEASUREMENTS OF SHORT PULSES 
OF DIRECT CURRENT. 
BY 


PHILIP F. ORDUNG.* + 


Yale University, 
1946. 


SYNOPSIS. 


Within recent years, engineering applications using pulses of direct current as short in 
duration as a microsecond and repeated several thousand times per second have become im- 
portant. Such pulses may have amplitudes that are greater than 100 amperes. This paper 
discusses the techniques and precautions which are essential for oscillographic observations of 
such pulses. Particular consideration is given to the problem of determination of the band- 
widths required for measurements of such pulses, of the construction and the calibration of 
resistors for metering such currents, and of properly terminating the coaxial transmission line 
used for connecting the metering element to the oscillograph. 


INTRODUCTION. 


With the advances and developments of recent years in the applica- 
tion of electron tubes has come an increasing need for accurate measure- 
ments of short pulses of direct current. The information which must 
be supplied by the measurements is the shape, amplitude, duration, and 
repetition rate of the pulse. Methods have been developed using meters 
for measuring the amplitude and the repetition rate of the pulse. How- 
ever, visual observation of the pulse by means of an oscillograph is the 
only method which can give accurate data concerning the shape and the 
duration of the pulse. The general arrangement of the components for 
the oscillographic method with which this paper is concerned is shown 
in Fig. 1. The pulse current, shown as ‘‘i” on the figure, flows through 
the metering element. The corresponding voltage pulse developed 
across the metering element is transmitted to and viewed on the oscil- 
lograph. For this type of measurement, the sweep of the oscillograph 
must be synchronized to the pulse being viewed and the interval of dura- 
tion of the sweep trace must be slightly greater than that of the pulse. 

The types of pulses that may be encountered in the applications of 
pulse techniques are of a great variety. Typical pulses of current have 
intervals of duration between a fraction of microsecond and several 
hundred microseconds and have amplitudes ranging from several am- 
peres to greater than one hundred amperes. One such pulse is shown 
in Fig. 5A. The shape of such pulses is approximately trapezoidal with 
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periods of rise and decay that are short in comparison to the pulse. 


duration. The repetition rates for these pulses range from several 
pulses per second to several thousand pulses per second. The problem 
of measuring such widely differing pulses is great, but the precautions 
which must be observed in their measurement are similar. 
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Fic. 1. Circuit for visual measurement of repeated pulses of current. 


BANDWIDTH FOR PULSE MEASUREMENTS. 


The measuring equipment in order to faithfully reproduce the shape 
of a repeated pulse must represent all of the frequency components 
contained in the pulse. For the purpose of determining the necessary 
bandwith for representing the pulse, Fourier Analysis is applied to the 
repeated square pulse shown in Fig. 2. J. is the period of the pulse, 
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Fic. 2. Wave of repeated square pulses of amplitude, A, 
period, 7x, and pulse width T>. 


and T> is the duration of the pulse. In functional form, the shape of 
one cycle of this pulse wave may be represented as follows: 


FO 
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where 
2rkt 


2 TR/2 

a= Cos —— dt. 2A 

é Tr f.. {2 ie ) - Tr ( 

The amplitudes of the frequencies required to represent the pulse are 

determined by substituting the functional form of the pulse, Eq. (1), 

into Eq. (2A). * For the case of the repeated pulse being discussed, the 
Fourier series sapeein ooe is 

a ee 2xrkt 

= — —C . 3 

f(t) Aunt T, + Usk Sin n rk = Cos T, (3) 

The angular frequency components contained in the repeated pulse 

are given in this series by the term 27k/Tp, where & is the order of the 

harmonic of the repetition frequency Fx. Using Eq. (3), the frequency 

spectrum shown in Fig. 3 is plotted. The envelope of the spectrum is 
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Fic. 3. Frequency spectrum corresponding to a repeated 
pulse wave in which Tz = 107p and Fp = 5Fp. 


shown by the dashed lines. Significant points in the envelope occur at 
frequencies that correspond to values of & such that 


UT r/Te = m  m=1,2,3,4,5 +--+ etc. (4) 


Since the term 27rk/Tp in Equation (3) expresses angular frequency, 
k/Tp is frequency. The particular value of this frequency, k/T., that 
is obtained when m = 1, is defined to be the pulse-width frequency Fp 
and has a period twice that of the pulse width, Fp = 1/2Tp. For the 
even integer values of m, the magnitude of the envelope is zero. For 
the odd harmonics of the pulse-width frequency (frequencies that corre- 
spond to the odd integer values of m), the amplitudes of the minor lobes 
of the spectral envelope are approximately maximum. Thus, changing 
the repetition rate of a repeated pulse has no effect on the shape of the 
envelope of the frequency spectrum ; however the number of frequencies 
represented in each lobe of the envelope is changed. 

For a square wave, Tp = T/2, the frequencies representing the 
pulse are odd harmonics of the pulse-width frequency. When the pulse- 
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width frequency is a multiple of the repetition frequency (Tr/Tp an 
integer), the spectrum includes the frequencies required to represent a 
square wave with the given pulse width and additional frequencies 
located symmetrically in each lobe of the spectral envelope about the 
square-wave frequencies. The additional frequencies eliminate the 
pulses in the square wave that are not contained in the wave having the 
reduced repetition rate. This suggests that the number of lobes in the 
envelope of the frequency spectrum required to approximate a pulse of 
given width to a desired accuracy is the same for a square wave and 
for repetition rates that are submultiples of the pulse-width frequency. 
Generalized, this argument extends to all repetition rates less than the 
pulse-width frequency provided that the pulse separation is such that 
wi(Tr — Tp) > 18. (See the Appendix.) 
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Fic. 4. Approximations to a square pulse of given width for two repetition rates— 
Fr = Fp, and Fr = Fp/2.85. The maximum frequencies included in the pass band for the 
approximations are as follows: (a) 2Fp; (b) 4Fp; (c) 6Fp; and (d) 8Fp. 


Approximations of a given square pulse for two repetition rates, one 
of which is that for a square wave, are shown in Fig. 4.. For both 
repetition rates, the inclusion of all frequencies to the eighth harmonic 
of the pulse-width frequency gives similar trapezoidal approximations 
having rise and decay periods less than 12 per cent. of the duration 
of the pulse and having less than 12 per cent. deviation in the tops of the 
pulses from the postulated shape. In the Appendix, the ratio of the 
rise time to the pulse duration is shown to be inversely proportional to 
the number of harmonics of the pulse-width frequency included in the 
reproducible band. Practically, all pulses that may be encountered 
have finite rise and decay periods which are seldom less than 5 per cent. 
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of the pulse width. Therefore, the bandwidth need seldom be greater 
tham the twentieth harmonic of the pulse-width frequency. 


CURRENT-METERING RESISTOR. 


For an oscillographic measurement of a current pulse, the shape of 
the voltage pulse developed across the metering resistor should be 
exactly similar to that of the current pulse. Therefore the metering 
resistor must have the same phase shift and amplitude characteristics 


Fic. 5. Photographs of a 6 ampere, 1 microsecond current pulse repeated 1000 pulses 
per second. 

(a) Pulse obtained with the corrected current-metering network connected with very short 
leads to the oscillograph. 

(b) Pulse obtained with the inductive current-metering network connected with very short 
leads to the oscillograph. 

(c) Pulse obtained with the corrected current-metering network connected with a 4-meter 
length of unterminated coaxial cable to the oscillograph. 

(d) Pulse obtained with the corrected current-metering network connected with a 4-meter 
length of unterminated coaxial cable to the oscillograph. 


for all frequencies contained in the band from zero to infinity. Such a 
resistor is physically impossible; but a network that approximates a 
perfect resistor over a\bandwidth sufficiently great to represent the 
current pulse shapes to within few per cent. error can be constructed. 
In the measurement of large pulses of current, several watts of average 
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power are dissipated in the metering network. For permanency of 
calibration and the dissipation of power, the noninductively wound type 
of wire resistor should be employed as the basic element in the metering 
network. The frequency characteristics of a typical commercially 
available noninductively wound 10-watt resistor are shown in Fig. 6B. 
The series residual inductance exhibited by the resistor is constant over 
the band of frequencies of immediate importance in pulse current 
measurements. The wave form of the voltage across this resistor when 
the trapezoidal pulse of current shown in Fig. 5A is passed through it is 
shown in Fig. 5B. The residual inductance is responsible for the oc- 
currence of the voltage spikes during the intervals that the pulse current 
has a high rate of change. 

The effects of the residual inductance are corrected in the network 
shown in Fig. 6A. When R: = R:, and Rik: = L;/C, the circuit is 
purely resistive with a resistance magnitude R;, that is independent of 
frequency. : The slight residual inductance in the resistance-capacitance 


y 


(b) 


Be 


Fic. 6A. Equivalent circuits of (a) the inductive metering resistor, and (b) the same resistor 
connected in a frequency correcting circuit. 


branch is the factor that limits the bandwidth for which the network 
approximates a pure resistance. The frequency characteristics for a 
typical network in which the inductive resistor previously discussed was 
used as the basic element are also shown in Fig. 6B. With this network 
the shape of current pulse shown in Fig. 5A was obtained. The exten- 
sion of the bandwidth for the network beyond that for the inductive 
resistor affected elimination of the voltage spike shown in Fig. 5B, but 
had no effect on the remainder of the pulse. Since the bandwidth ex- 
tends to 30 megacycles, this particular network may be used to view 
pulses of current as short as } microsecond with a rise time that is 8 
per cent. of the pulse width. 

Practical application of this type of network is made feasible by the 
relative power dissipation requirements of the inductive and the capaci- 
tive branches under pulse conditions. During the pulse rise the con- 
denser in the network charges with an efficiency of 50 per cent. through 
its series resistor R:, and during the pulse decay discharges through this 
resistor. The energy dissipated by R. during the pulse rise is equal to 
the maximum energy that is stored in the condenser. An equal quantity: 
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of energy is dissipated in R. during the pulse decay. Thus the total 
energy dissipated in Rz per pulse is twice the maximum energy stored 
in the condenser. The average power loss in R: is the product of the 
repetition rate and twice the maximum energy stored in the condenser 
during the pulse. If the network (b) shown in Fig. 6A were used to view 
a 20-ampere 1-microsecond square pulse having a repetition rate of 
1000 pulses per second, the power loss in R, would be (CE*) X (Repeti- 


—— INDUCTIVE RESISTOR 
——-CORRECTED RESISTOR 


0 


lO IS 
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Fic. 6B. Equivalent series-circuit frequency characteristics of the inductive 
metering resistor and of the corrected resistor. 


tion Rate) or 0.122 watt. The remainder of the power input to the net- 
work, 5.87 watts, would be dissipated in the inductive resistor, R;. For 
such a low power loss Re may be a 1-watt composition resistor having a 
sufficiently small residual inductance that the R. — C branch of the 
network may effectively compensate for the inductance in R. 


TRANSMISSION OF THE PULSE FROM THE RESISTOR TO THE OSCILLOGRAPH. 


From the standpoint of measurement technique the leads carrying 
the pulse current to the metering resistor should be short, and one 
terminal of the metering resistor should be as near to ground potential 
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as possible. Therefore lengthy connections from the metering resistor 
to the deflection plates of the oscillograph are usually necessary. To 
keep the viewing circuit as free as possible from extraneous effects, the 
connections should be made with coaxial cable. When measurements 
are to be made on short pulses, considerable care must be observed in the 
load termination of the coaxial cable. The efficacy of the termination 
of a given cable may be evaluated from measurements of its input 
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Fic. 7. Frequency characteristics of a 4-meter length of a 90-ohm coaxial cable with 
different terminations—(a) unterminated; (b) terminated with a single resistance; (c) termi- 
nated with the type of network shown in Figure 8. 


impedance. When properly terminated, the input impedance is the 
characteristic impedance. When improperly terminated, the cable 
causes transient reflections which are superimposed upon the observed 
pulse and tend to obscure its shape. The amplitude of the reflections 
are indicated by and may be evaluated from the deviation in phase 
and magnitude of the input impedance from the characteristic im- 
pedance of the cable. The observed pulse when a 4-meter length of 
coaxial cable without termination (see Fig. 7 for the input impedance of 
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the cable) is used in the measurement of a trapezoidal current pulse as 
shown in Fig. 5C. The transient reflections make the true shape of the 
one microsecond pulse completely indistinguishable. 

When the coaxial cable is properly terminated in its own character- 
istic impedance, reflections of this type are completely eliminated. In 
coaxial cables suited for very high frequency transmission, the attenua- 
tion at the frequencies used in pulse measurements is sufficiently low 
that the characteristic impedance is nearly constant in magnitude and 
practically resistive in nature. Therefore the termination should be 
purely resistive. However, the input capacitance of the oscillograph, 
Cy = 19 wyf in this case, is part of the termination. For the frequencies 
needed to represent pulses one microsecond or greater in duration, the 
shunting effects of the input capacitance are sufficiently small so that a 
simple resistance termination may be used (see Fig. 7). For shorter 
pulse lengths a network which compensates for the input capacitance, 
Co, of the oscillograph is necessary. A typical network electrically simi- 
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Fic. 8. Coaxial cable termination network for pulse lengths that the input 
capacitance, Co, of the oscillograph is important. 


lar to the type used in the current metering resistor is shown in Fig. 8. 
When R; = Rz and R:R, = L/Co, the impedance of the network is 
purely resistive and of magnitude R;. Since the power dissipated in 
this network is very small, low wattage components that differ little 
from the ideal over the pulse bandwidths may be used. The improve- 
ment in the transmission bandwidth for the pulse cable when this net- 
work is used for termination is shown by the solid line in Fig. 7. Since 
the bandwidth extends to 30 megacycles, the cable may be used in the 
measurement of pulses are short as } microsecond. 


SYSTEM CALIBRATION. 


For measurement of the amplitude of the pulse, the deflection sensi- 
tivity of the oscillograph and the equivalent resistance of the metering 
network and terminated coaxial cable are required. The deflection 
sensitivity is measured directly atthe oscillograph by observing the 
deflection of the beam for a known d-c. voltage. The equivalent re- 
sistance of the network and cable that is presented to the pulse is 
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measured at the pulse-current terminals of the network. This value 
measured with a Wheatstone Bridge is accurate for the bandwidth of the 
metering system and is equivalent to the resistance of the metering 
network and characteristic impedance of the terminated cable consid- 
ered in parallel. When short lengths of coaxial cable are used, the 
attenuation of the pulse between the metering network and the oscil- 
lograph is negligible. Therefore the current sensitivity of the metering 
system in amperes per unit of cathode-ray deflection is the quotient of 
the deflection sensitivity of the oscillograph and the equivalent resist- 
ance of the metering system. 

Measurement of the duration and specification of the shape of a 
pulse requires that the sweep of the oscillograph be calibrated against 
time. A circuit for calibration is shown in Fig. 9. The synchronizing 
pulse-generator starts the pulse generator and triggers the sweep of the 
oscillograph. The output of the calibration pulse generator is a series of 
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Fic. 9. Circuit for calibration of the sweep time of the oscillograph. 


pulses that appear on the screen of the oscillograph and that are short 
in duration compared to the sweep. The time interval between the 
pulses shown on the screen of the oscillograph is the period of the funda- 
mental frequency of the calibration pulse generator. If the repetition 
frequency of the synchronizing pulse generator is low, less than several 
hundred impulses per second, and if the interval of operation of the 
calibration pulse generator is greater than half of the interval between 
synchronizing impulses from the synchronizing pulse generator, the 
frequency of the calibration pulse generator may be measured in a 
conventional manner with a calibrated radio receiver. Thus the sweep 
of the oscillograph may be calibrated against time. 


APPENDIX. 


The bandwidth required for a pulse to be represented with a given 
accuracy may be obtained more rigorously by the use of Fourier Trans- 
forms. In this derivation the fundamental Fourier Transform rela- 
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tionships ! are assumed as follows: 


g(w) = 5- f* sera, 


Hi) = [* g(w)er'de, 


The time function f(¢) defines the disturbance that is to be analyzed. 
g(w) is the corresponding frequency function that reproduces the dis- 
turbance. In this derivation f(t), shown in Fig. 10A, is a repeated 


I 


Geo , 


Fic. 10A. Repeated pulse wave f(#) initiated at ¢ = 0 and continuing thereafter. 
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Fic. 10B. Response function F(t) showing the first two pulses of f(t) for the particular case 
that Tr = 20, Tp = 10, and the bandwith, w; = x. 
pulse wave initiated at ¢ = 0. In the mathematical form which is 
used in the transforms, the time function is defined as follows: 
f(t) = 0, t< 0, rTr+Tp<t< (r+ 1)Tp, 


7 
f(t) = Ae-**, rI'p <t<rTr+Tp, r=0,1,2,3-:-. ) 
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The integrand in the time integral, Eq. (5), must be convergent in order 
for the integral to exist. The factor e-** was introduced into the ampli- 
tude function of the wave to satisfy the convergence requirement of the 
time integral. Since the constant a may be chosen small, the effect 
of the e~«‘ term on the first few pulses in the repeated wave is negligible. 
The frequency function g(w) that reproduces this wave is obtained when 
the function f(t) defined in Eq. (7) is substituted into the infinite time 
integral as follows: 
1 Tp Bee ae rTR+TP 
g(w) tie ~f Ae ati») td} + leas > f Ae~‘«ti+)td} 
2 0 ar rTR 


us r=] 
ie 4] 1 Ba e~(et+io Tp 
~ Qrla t+ jw a + jw 


+ G(w), (8) 


where 
rTrR+Tp 


G(w) = a > f Ae~‘«+i~) dt, 
2r ray VTTR 

Here the frequency function is written as the sum of two parts, that 

which arises from integration over the first pulse and that which arises 

from integration over the remainder of the pulses. This form of the 

frequency function is retained throughout the derivation. 

If this frequency function were transmitted through a perfect 
network, one which has either zero or linear phase shift and unity 
transmission for all frequencies, the time function would be reproduced 
but with the exception of a Gibbs effect at the extremities of each pulse. 
For a practical network, the pass band is finite and may be defined as 
—w: <w<:. In the pass band the network is presumed to have 
ideal characteristics—linear phase shift and unity transmission. Out- 
side of the pass band the network has zero transmission. In practical 
networks used for pulse measurements the frequency characteristics 
are sufficiently similar to the idealized characteristics that calculations 
based upon the ideal are useful for obtaining the response of an actual 
network to a pulse. 

The response F(t) of the network to the frequency function g(w) 
when the pass band is limited to —w; < w < w: is determined by sub- 
stituting the product of g(w) and the frequency characteristics of the 
network into the frequency integral, Eq. (6), with the limits of integra- 
tion restricted to the pass band. Thus 


F(t) == ss g(w)ete\'-dy 


1 e-(atiw) Tp 


atjo a + jw 


eie(t—todg) 


+f" G(w)ei#-dw, (9) 
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The term, e~‘*'e, where ¢ is a constant, that appears in the integral 
expresses the linear phase-shift and unity-response characteristics for 
the frequencies in the pass band of the network. By a substitution of 
t — t for t in Eq. (6), the linear phase shift term e-‘** is shown to have 
the effect of delaying the network response function F(t) by a time éo. 
Since the magnitude of the time delay is of no importance in the re- 
mainder of the derivation, ¢) is not evaluated. In Eq. (9), the integrand 
for negative frequencies is conjugate to the integrand for corresponding 
positive frequencies. Thus by writing the integral in Eq. (9) in the 


@1 0 @1 
form f = f + f and by substituting —w for w in the integral for 
—@1 —@) 


the negative frequencies, the response of the network may be written 
in the form 


as #1 4 £( oe) oO 
Fé) A gg ER) erin 


f° —(a-—jw)Tp 
A A(t if - ) e—i(t—todgy 
ai jo a— 7JW 


+ | Glw)ei#'-dw, (10) 


—wl 


Eliminating the complex functions in the integrals gives 


A a 
F(t) = — To [Cos w(t — to) 


0 a 
— e~*TP Cos w{(t — to) — Tp} |dw 
+i [" eset — to) 


a? + w? 


dw 


— to) mt T p}dw 


+ | Glw)ei-de, (11) 


—w) 


The damping term e~*' was originally introduced to provide the needed 
convergence for the time integral. For small values of a, the first few 
pulses of the original time function f(#) are not appreciably affected. 
As @ is decreased, the number of pulses unappreciably effected by the 
e~* term becomes larger. As the Lim a — 0, f(#) approaches a square 
pulse wave beginning at ¢ = 0 and continuing for an infinite time. In 
the remainder of the analysis, the limiting condition of a is used. 
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As Lim a — 0, the following conditions apply to the terms in Eq. (11): 


Lim e: at rs [Cos w(t—t)) —e-@TP Cos w{(t—to) —Tp} Jdw=0 
ao Y9 @?+a’ 
tim | —*—Sin ott tide = Si 6,02) 


an J0 at+u? 


; vas ae ‘ ; r 
Lim , Gat nin wl (t-te) — Tr}dw=Si wil (tt) — Tr}. 


Thus the response of the network becomes 


ft) = 4 [Si w(t — 4) — Siwaf(t — &) — Tr}] 


+o1 
+ Lim f G(w)ei#'-dw, (13) 
a—0 —@i 

The functions,? 1/2 Si wi(t — to), 1/4 Si wi{(t — to.) — Tp}, and the 
first two pulses represented by Eq. (13) are shown plotted in Fig. 10B. 
For values of |wi(¢ — to)| > 18, the value of 1/x Si wi(t — to) ~ +}. 
In repeated pulse waves where the time of separation between pulses is 
such that w:(Tr — Tp) > 18, each individual pulse may be described 
with negligible interference from the terms for other pulses by a term 
similar in form to A/a[Si wi(t — to.) — Siwi{(¢ — to.) — Tp}. Under 
these conditions, the first term in Eq. (13) describes the first pulse in 
F(t). The remainder of the pulses in F(t) are similar in form but de- 
layed in time from the first and are contained in the integral term of the 
equation. Thus the accuracy of the approximation obtained by limit- 
ing the bandwidth to —w: < w < w: may be determined from the 
behavior of the first term in Eq. (13). For bandwidths such that the 
shape of the pulse is approximately trapezoidal, the slope of the tangent 
to the pulse at ¢ — t) = 0 is about average for the rising edge of the 
pulse, and has the equation 
< F(t) 5 [Si wi(t—to) —Si wif (t—to) — Tp} ] Ph 
The time of rise, tise, of the pulse is calculated from the slope of the 
tangent to the pulse at ¢ — ¢t, = 0 as follows: 


hin = 5 . — (15) 
. W1 
dt F(t)) suo 


{ 


(14) 


t—te =0 T 


However, the accuracy of a given approximation to a pulse should be 
expressed in terms of the ratio K of the rise time to the duration of the 
pulse. 
brise T 1 
=O = . 16) 
<“t, ote Wee \ 
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In order to obtain a convenient expression for the bandwidth needed to 
give a desired approximation to a pulse, F; is replaced nFp, where n 
is the harmonic of the pulse-width frequency Fp». With this substitu- 
tion, Eq. (16) reduces to 


Xs 


Therefore the ratio of the time of rise to the pulse width for a given re- 
peated pulse wave is inversely proportional to the number of harmonics 
of the pulse-width frequency that are included in the pass band, provided 
that the separation between pulses is such that w:(Tr — Tp) > 18. 
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LIST OF SYMBOLS. 


a, Amplitude coefficient of the kth harmonic of Fr. 
A Amplitude of the pulse in the repeated pulse wave. 
C)  Oscillograph input capacitance. 
f(t) The repeated pulse wave, defined differently in body and Ap- 
pendix. 
F,; Frequency corresponding to w. 
Fp  Pulse-width frequency. Has period equal to 2T >. 
Fe Basic pulse repetition frequency. Has period equal to Tx. 
F(t) The response function of the network when the bandwidth is 
limited to w,. 
g(w) The frequency function corresponding to the f(t) used in the 
transforms. 
G(w) A portion of g(w) in the Appendix. 
k Number of the harmonic of Fe. 
Ratio of the rise time to the pulse width. 
Series of positive integers which are involved in envelope of 
spectrum. 
Number of harmonics of pulse-width frequency included in pass 
band. 
r The number of pulses that has occurred counted from the first. 
R,, Ro, L, C The constant in the measuring network. 
t Time. 
ty Time delay for a pulse. 
T» Duration of the pulse. 
Tr Interval of time between successive pulses in the pulse wave. 
» Angular frequency. 
», Limit of the pass band for the measuring networks. 
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Frozen Juice Like Fresh.—One of the very best of the orange juice products 
is a semi-concentrated juice stored at zero temperature and made ready for 
drinking by diluting with three parts of water to one of juice. 

This is the opinion of workers at the U.S. Department of Agriculture's 
Winter Haven, Florida, Laboratory, who haveh ad opportunity to taste the 
products of many systems of processing. It was developed at the Laboratory 
in cooperation with the Florida Citrus Commission. 

“In this method,’”’ says M. K. Veldhuis of the Bureau of Agricultural and 
Industrial Chemistry, ‘orange juice is vacuum-concentrated to about 5-fold. 
Then fresh juice is added to reduce the mixture to a 4-fold concentration. 

“The addition of fresh juices,’ says Veldhuis, ‘‘replaces much of the flavor 
lost in the vacuum concentration process.” 

This 4-fold concentrate is then quick-frozen and marketed through frozen 
food channels. ‘‘At the usual zero storage temperature, the product is little 
harder than a slush and can be mixed readily with three volumes of tap water 
to quickly prepare a cold drink hardly distinguishable from fresh juice but 
with a decided gain from reducing costs for transportation and storage.” 


R. H. O. 


X-Raying Oil Fields. (Electrical Engineering, Vol. 66, No. 7.)—A method 
of X-raying oil fields by remote operation from its research laboratories has 
been developed by Gulf Oil Corporation. The new development helps solve 
the riddle of how oil flows through the sand and rocks under any particular 
field, and how it is affected by underground water and gas. 

Since X-raying entire oil fields at their sites would be impractical, the fields 
were brought to the laboratory. Small 3-inch by 1-inch samples, or cores, from 
the drill holes are made to reproduce in miniature the layer of rock or sand from 
which they were taken. Reactions studied within the samples give a picture of 
flow conditions in a stratum perhaps hundreds of feet thick and miles in extent. 

To make such studies, the core is subjected to artificial pressures and sat- 
urations of oil, gas, and water. The progress of these elements in the core is 
gauged by an X-ray device for determining permeability-saturation. The X- 
ray beams follow the reaction by means of an opaque tracer mixed with the 
liquid or gas. 

First step of the researcher is to establish a basis for comparison by studies 
of the core at 100 per cent and at zero saturation. Then an extensive series of 
flow experiments, reproducing conditions which might be created by various 
recovery methods, must be made. 

Effects of gravity, capillary attraction, and the amounts of oil, gas, and 
water already in the sand are charted. These procedures must be repeated 
sometimes on as many as 10 or 11 cores taken from the different layers through 
which the well extends. 

By correlating such data, the laboratory can determine how natural press- 
ure, or artificial gas injection, or water flooding will move oil through various 
strata of the field being studied toward well shafts. Operating methods can be 
planned accordingly to assure the greatest yield, and the extent of the yield 
forecast. 

Permeability-saturation determination by X-rays represents an extension 
of the electrolytic model oil field technique developed for predicting generalized 
oil field behavior and guiding placement of wells. R. HO 
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RIGOROUS AND APPROXIMATE TREATMENT OF LONG LINE 
TRANSMISSION PROBLEMS BY HYPERBOLIC FUNCTIONS.* 


BY 
S. GERSZONOWICZ. 


Professor of Electrical Engineering, University of Montevideo. 
I. INTRODUCTION. 


It is generally agreed that in rigorous treatment of long transmission 
lines hyperbolic functions must be employed. In order to facilitate the 
application of these functions, charts and tables have been prepared; 
the last chart, due to Woodruff, covers the values of 6? = ZY between 
0.05 and 0.36, corresponding to the lengths of aerial lines between ap- 
proximately 100 and 300 miles at 60 c/s, because 6/L is practically con- 
stant and equal to 0.0018-0.0022 per mile at 60 c/s. 

Almost all authors state that when using the expansions in series of 
hyperbolic functions it is possible to obtain any precision desired by 
taking a sufficient number of terms. They consider however the mathe- 
matical precision only; many of them (for example, Nesbitt, Wagner 
and Evans, etc.) give in their books examples of computation of suc- 
cessive terms in order to show that these terms diminish very rapidly, 
so that quite soon they do not affect the first, say four, significant figures. 

A reference—rather vague—to the physical problem is found in G. 
D. McCann’s excellent chapter No. 4 in the ‘Westinghouse Trans- 
mission aad Distribution Reference Book”’ (1944). McCann fixes at 
0.5 per cent. the mathematical precision of the developments, because 
“the resistance, inductance and capacitance of a line can rarely be known 
within 3 or 4 per cent. and probably never within 1 per cent. This is 
due to conductor sag, its variation with differents spans and the varia- 
tion that exists in conductor spacing, together with the effect of tem- 
perature upon conductor resistivity and sag.”’ 

In this paper we shall treat the physical aspect of the problem more 
comprehensively. We must consider two points: 


(I) There is no physical significance in adding a new term of a 
development unless this term is large enough as compared to the maxi- 
mum error committed in the value of the sum of preceding terms. We 
shall admit that it is so if a consecutive term is larger than 0.25 of the 
maximum error in the sum of the preceding terms. Of course, a value 
other than 0.25 may be chosen without affecting the generality of the 
treatment; 

* Lecture given at the Graduate School of Engineering, Harvard University, on April 10th, 
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(Il) Even if it makes sense to add a new term in compliance with 
the consideration (1), its addition may not be worthwhile in view of the 
limited accuracy within which a given magnitude has to be known. 


Asa corollary of (II) we shall examine the problem under which con- 
ditions a line may be considered short, that is of negligible capacitance. 
The textbooks agree that an overhead line may be considered short 
when it is 20 or 30 miles long; some say that the effects of capacitance 
may be neglected whatever the length of line, if the voltage is lower 
than 40 kV. We shall see that the problem cannot be stated in this 
simple manner. 

It is easy to see, applying to the term in @* the same treatment that 
we shall use in respect to the terms in 6? and @', that the term in 6° may 
be omitted for present usual lengths of lines. Consequently we shall 
limit all developments to the terms in @ and below. 


II. ERRORS ON THE CONSTANTS OF A LINE. 
Let, per mile, 
r—the resistance, in ohms, 
lw—the inductive reactance, in ohms, 
cw—the capacitive susceptance, in mhos. 


We shall neglect the conductance g, ordinarily very small. We have 
Z = Aly = Vr’ + PotL/arctan lo/r, 
Y = Y/x/2 = cwL/90°, 


where L is the length of the line in miles. 
The maximum errors on Z, y and Y are 


+ Ae) i sin’ y a, 
A(lw) 


Ie 


+ - x) sin y COs Y, 


and are due to the variable atmospheric conditions, the variation o! 


frequency and of the load, and the approximation of the formulas used 
for the calculation of the constants, or the errors in the measurements 


of these constants. 
Il. THE EQUIVALENT ~ CIRCUIT. 
We have 
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A. The Development of E 


We have, limiting the development to the term in @, 


Z.= 2(1+2+4) (6) 


and the projections of Z., M on Z and N on the normal to Z, are (Fig. 1): 


‘ OF 6 
uw =z(1-"siny— cos2y), 


e @ 
tany = N/M= 6 COS ¥ + 75, 


The modulus and argument of Z, are 


peeing =~ M(1 + $ tan’ y) 
cosy 


j -.. 6 |, 
Pe ae + ( 
= z| 1 6 SiN Y — 799 ( cos 2y — $ cos’ y) (9) 


4 


180 sin 27. (10) 


@2 
+-—cosy + 
6 
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If the use of the second term, (Z@ sin y)/6, in (9) has to have any Cor 
physical meaning, we must have, according to the consideration (1): 


: .. * 
tAZ < (Z@ sin y)/6; hence 6 > —- a (11) @@ but als 
2siny Z and ou 


accoun 


In examining the use of the third term, we should compare it with the 
maximum error in Z — (Z@ sin y)/6, but in an analysis of errors it will 

be sufficiently correct to compare it with the maximum error in Z only, We 
(Z@ sin y)/6 being relatively small with respect to Z. _Consequently, 

the use of the term in @ in (9) has no physical meaning unless 


SOE Sa. AES. 2 
— |cos 2y — $ cos? y| Z 


(12 


Let us consider now the argument y’. The use of the second term in 
(10) has a meaning only if 


3 
> — Ay (13 
2 cos 


and the use of the term in @ in (10) has a meaning if 


(14 


To see the actual significance of the relations (11)—(14), and of the 
similar ones which we are going to establish for other magnitudes, we 
have to assume numerical values for errors. We shall admit 


A(cw) _ A(lw) _ AL 


= 0.01 
cw lw L 


and Ar/r = 0.06, what corresponds to a 15° C. variation in tempera- 
ture. The relations (11 and 13) become: 


> 30.06 cos* y + 0.01 sin* y + 0.01 


: (Ila 
2 sin ¥ 


> 


# > 0.105 sin y. (13a 


The relations (11a) and (13a) show that the lengths L’ (modulus) and 
L"” (argument) for which the terms in 6? have to be considered in (9) 
and (10), respectively, depend on y. Making @ = 2.10-*L, we find The sec 


modulus (lla): L’ = 85 miles for y = 90°, and 155 miles for 7 = 45°; 
argument (13a): L’’ = 160 miles for y = 90°, and 135 miles for y = 45°. 


Similarly, the relation (12) gives and the 
L’ = 450 miles for y = 90°, and 600 miles for y = 45°. 


The relation (14) gives L’’ = 550 miles independently of y. 
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Consequently for present lines the term in @ should not be taken into 


TP v 6 
account. It not only means that it is sufficient to make Z, = 2( 1+ *). 


but also that, when calculating Z., we do not have to correct for 1/cos y, 
and ought to write. Z, = M. 
‘h the 
it will 
only, 


ently, Y #2 4 
Y, = A i+ D sin y — ~~~ (cos 2y — +s cos’ y) |. (15) 


B. The Development of y.. 


We find, in similar manner as for Z, (Fig. 2): 


2 120 


— 


? 


_ cosy 7... Fe 
tan M7 = - a0 ¥ ( 1 + 60 = nee , (16) 


pera- 


(lla 


) and Fic. 2 


The second term in (15) will be used if 


and the third one if 


> - 
~ |cos 2y — 5 cos? y 
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Let us consider now the argument of F,., (90° — ¥). In view that we 
assumed g = 0, the term in 6 in (16) must always be considered, as 
far as only the consideration (I) is concerned. 

The second term in (16) has to be considered only if 


ay. Bf ae 
2 alee 


C. 


Until now we have considered Z, and Y, alone, but generally it is 
necessary to consider them in combination with themselves and with 
other magnitudes. The conclusions obtained before may not be valid, 
because of possibility of compensation and of the presence of new terms. 

We are interested in obtaining correct values of the terms in 6 and 
6* in the developments corresponding to various magnitudes. Note 


: 5. 


=2 120 NNIVAA 


that for expression of a form Z,*¥.’, where a and b are positive in- 


6 ree e 
tegers or zero, it will be sufficient to make Z, = Z ( 1 + *) and 


Y Ce — hay ae 
y, =  ( 1 — 5) ; the coefficient of the term in & will be always 


exact, and the one of the term in @ will be exact if a > 0 and b > 0; 
if a = 0 or b = O, this coefficient will be too small in (a + 6)/120. 

If we make Z, = Z and Y, = JY, the coefficient of & will be exact 
only if a > O and b > 0. 


IV. LINE AT NO LOAD. RELATION OF THE VOLTAGES. 


We have (Fig. 3): 


U so -. - 
—=1+4+2Z.7,=A = A/é, (20) 


Uro fe 


where A is one of the constants A, B, C and D = A of the equivalent 
four terminal network. 
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Let us make in (20) Z, 


we obtain 


where the coefficient of # is exact, according to what we said above. 
The modulus is  ° 


Us -.. . 
at tee tie l1—-~sny -= g \c0s 27 — 3 cos? > y) (22) 


Uro 2 


and there is a meaning in using the term in @* if this term is larger than 
0.25 of the error in the term in 6°, that is, if 


3 sin y [te 4. 8 ey AZ. (23) 


2> 
, la Y 


—~ 1 + cos? y 


;A l af 
: “a += — os — = 0.04, for y = 90° the length must be greater 
than veg miles, saat a y = 45° greater than 120 miles. 

The value of A is obtained with great accuracy by means of (22); 
for example, if Z = 120 miles, we may neglect the term in 6* (which 
means here to make Z, = Z, Y, = Y/2, and cos 6 = 1) and the maxi- 
mum error on A will be between ¢ and 1.25 ¢«, with 

@ si A(l ab . oY 
i sn (lw) vt (24) 


2 lw bes elit 


For L = .120 miles, and aoe) + — =) 2 —- = ().04,e 0.001. But 


normally it is sufficient that A have a pte less accuracy. Indeed, if 
we determine A by measuring Uso and Uro on the respective switch- 
board voltmeters, it is impossible to expect to obtain A within a maxi- 
mum error less than 2 per cent. In such a case it is useless to take into 
account the term (6? sin y)/2 if it is smaller than 0.02, that is, if sin y = 1, 
for lines of less than 100 miles. Then, under the conditions just exposed 
aline up to 100 miles has to be considered as “‘short.”’ 


V. LINE AT NO LOAD. CURRENT AND POWER AT THE SENDING END. 
A. Current at the Sending End. 


We have (Fig. 3) 
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Let us make 2, = 2(1+5) nd P= > (1- ‘ 
S ake 42. >= 6 a oe: = D) 12 


-7(1-74+¢ ) (26) 


with an error in the coefficient of 6* equal to 1/120; the exact value oj 
this coefficient is 2/15. 
The modulus of Is0/U so is . 


a v| | + 5 sin y — (cos 2y — fa cost y) | 


and there only is a meaning in using the term in 6 if 
15 1 | 
<a ae 
—~ 8 |cos 2y — ¥s cos? y| Y 


that is, if AY/ Y = 0.02, 220 miles for y = 90° and 320 miles for 
For any length shorter than defined by equation (28) we may write 


Iso ‘eg v ( 1 +" sin y ) 
Uso Pare) 


with a maximum error between AY/Y and 1.25 AY/Y. 

To neglect the term in @ is not equivalent here to make Z, = Z ani 
Y.= Y/2. If we limit the developments of Z, and FY, to their first 
terms, we find 


(30 
with a supplementary negative error (@ sin y)/12. 


If it is sufficient to obtain J s9/U s» with an even less accuracy, we 
may write, neglecting Z.¥. compared to 1, 


Iso > v( . ) 7 
gee Pe aoe 3] 
U ec 2¥ 1 + = sin y (3 


with a supplementary negative error (@ sin y)/4, or 


I so ‘ 
7 Be } 
U So 


with a supplementary negative error of (@ sin y)/3. We have now to 
investigate when a line may be considered as “short” from the point 


of view of the no-load current. In order to give a physical significance 


to this problem we must refer J so to the load current J,. Assume that 
the line is transmitting its ‘‘natural’’ power, magnitude introduced by 


Riidenberg to characterize a load for which the line operates at nearly 
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its maximum efficiency. © Under these conditions 


Z sin 
Ps 303 / —— 
and hence 


Fas 
lue of I, = Urcos oad — 


(26) 


where COS gp is the power factor of the receiver. 
From (32) and (34) 


Tso _ , o> Y 2.10-3L 
I, Gr 


Z and 


r first 


(30 


a 


Fic. 4. 


with Lin miles. UsoVsin y/Uzx cos ¢p is nearly equal to 1; consequently 
we find that to call a line of a length less than 20 miles ‘‘short’”’ means to 
neglect a no-load current equal approximately to 5 per cent of the full- 


w load current.. 
po in B. Power at the sending end. 


yw to 


sance : ohne . Oe Di 
that We have, calling y the angle between J s0/U so and Y (Fig. 4): 

d by Pso _ Iso. 6° a: | 
early Uso? = 5 ge siny = Y ( 3 cos y + 15 sin 2v) (36) 
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and the use of the term in @ has a meaning if 


he -(* aL, 24Y = 
=“16siny\r '° L y /’ sig 


that is, if Ar/r + AL/L + 2AY/Y = 0.11 and sin y = 1, for lengths 
greater than 90 miles. 


VI. LINE IN SHORT-CIRCUIT. CURRENT AT THE SENDING END. 
We have 


(38) 


= 2 7 zs F & : 
Let us make Z, = Z ( ‘ = — } ; we obtain 


Iss ( j f 
Iss _1 aR 39) 
Pee | Cae aa is 


with an error of 6*/120; the exact value of the coefficient of # is — 1/45. 
The modulus is 
ten 8 r . 6 
—— «= —| | ——g@n 7 + = (cos2 $ cos* 
ia Be) eat een 
and the use of the term in @ only has a meaning if 
AZ 
ar eee pee (41 
+ cos’ vy Z 


that is, if Ar/r = 0.06, A(Jw)/lw = 0.01 and AL/L-= 0.01, if the length 
is greater than 350 miles for y = 90° and 400 miles for y = 45°. We 
may then write, with an accuracy between AZ/Z and 1.25 AZ/Z, 


fea 9 ee oe 
iss = 5 (1-siny ). (42 


If we may tolerate a larger error, we can make Z, = Zand Y, = Y/2, 
obtaining 
i r : 
LF Saees« (43) 
7( 2 98 y) , 


the supplementary error being negative and equal to (@ sin +) /6. 
We may also simplify (38) neglecting Z,.¥. compared to 1; we find 


1 1 (  . ) (44) 
in ea ee (44) 
7 Z . 6 re 


with a supplementary positive error (@ sin y)/2. 


of the ; 


Jan. 194 


Iss/U. 
A(lw) /1 
than al 
be 


where 


Let 


For cos 
we ma’ 


For 


express 


Jan., 1948.] Lone LINE TRANSMISSION PROBLEMS. 63 


Finally, let us observe that the simplest expression for Iss/U ss, 


t46 .. 1 

Cn ee he 
is equivalent to the consideration of the line as “short.” It leads 
to a supplementary positive error of (@ sin y)/3. If we want to find 
Iss/Uss within a total accuracy of 10 per cent., with Ar/r = 0.06, 
A(lw)/lw = 0.01, AL/L = 0.01 and sin y = 1, the length should be less 
than about 240 miles. For 300 miles, using (45), the total error would 
be = 14 per cent. if y = 90°, 13 per cent. of y = 45°. 


VII. LINE AT LOAD. 


A. Relation of the Voltages. 


We have 


_ wit SP 4 27. w 2 + &f. (46) 
R 


where Y, is the admittance of the load. 
8 es *( ) ie 
Let us make 2, = 2( 1+ and 7+ iach ; we find 


Ory (14%) 414%) 
ty NET PERE FT — 
If the power transmitted is the natural power of the line, we have 


y : - cos OR 1 
Y,. = \ COS ¢r, and — = ——-: 48) 
Z sin y a Y vsin y 8 ( 


For cos gr & 1 = sin y, and a length of 300 miles, Y./Y = 1.6, so that 
we may admit, only for the purpose of discussion of errors, 


L oe 8, 
tg Met g My yi #2 


For the same purpose we may neglect the correction of the cosine 
of the angle between 1 and Us/Ur, and find the following approximate 
expression for the modulus of Us/Ur: 

Us 


Z oe . 
so a 1 + ZY, Cos (y — gn) — —— ain (27 — 92), (49) 
Ur , 


and the use of the third term in (49) only has a meaning if 
AL 


3 
fe > shai 
~ 2 sin (2y — —| L oar = oe) 


+ = COs Y COS gr + ae) sin 7 SiN gp | » (50) 
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for example, if y = 70°, gz & 25°, Ar/r = 0.06, A(lw)/dw = 0.01 = AL/L, 
for L > 110 miles. For lengths less than the one defined by (50), 
we may make Z, = Z and Y, = Y/2 in (46), which is equivalent to 
neglecting the third term in (47). For greater "egeths all depends on 
the accuracy within which we want to know the ty ©'ation U s/Ur. For 
a line of 300 miles, to neglect the term in @ in (~"4 is equivalent to a 
positive error equal to 0.02 Y./Y, that is, about. es cent. if the natural 
power is transmitted. It must be observed th* bis error is propor- 
tional to the fourth power of the length, L‘. 


B. Relation of the Currents. 
NV e 
We hav 


= =1+2.7. + 2.7.(2 + 2.7.) = A+ C2Z., (51) 


2) +5(a+a) 


If the line transmits the natural power 


x Vsin 71 
COS gp 0’ 


that is, for a line of 300 miles, siny 2&1 and cos gr = 1, Z,/Z = 
1, 6/g2 — vy, so that we may still make for the purpose of the discussion 


of errors 


oe Ze Ze 
pte ea tee 
Furthermore we may nee for the same purpose the correction of the 


cosine of the angle between 1 and Js/Ix, and find the following approxi- 
mate expression for the modulus of Is/Ir: 


s : ZY? si 
Is =1—Z,Ysin gr — wc al cos (y + ¢r). (54) 
Ir 6 


Here the third term has a meaning if 
sin QR 
Y |cos (y + ¢a)|’ 


that is, if AY/Y = 0.02, y = 70° and ger & 25°, for lengths greater than 
200 miles. Furthermore it may often be neglected because we only 
need to know J s/Jz within a limited accuracy. For a line of 300 miles, 
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and for cos (y + ¢r) as high as 0.5 the error in neglecting the third 
term, thus making Z, = Z and Y, = Y/2 in (51), is positive and lower 
than 2 per cent. 

VIII. LINE AT LOAD.! POWER CIRCLES AT SENDING AND RECEIVING ENDS. 


We have 
(56) 
and 
(57) 


The modulus of the re'!ation A/B has been studied in VI. We found 
that in practically all cases it is sufficient to take 


A 1 e . 1 lwL VY 
£=5(1-Fsnr) =5(1- 3 ): (58) 


The argument of A/B is 
@? a 
vy" =-—yt+ 3 COS y(1 + 756 sin y), (59) 


where the use of the term in 6* only has a meaning if 


45] A(lw) =a 
ae Gi, “ph ie AEN Ua 
e2 | lw age ies 


that is, for A(lw)/lw + Ar/r = 0.07, for L > 280 miles. 
We may then write for all usual cases 


6? ay. 
y= y+ 50087 = — y+ —- (61) 


: mae fe | Ur 1 ; ; 
For the radii, UeZ, and UZ.’ according to the study done in II], 


is practically always sufficient to write 


mer). (62) 


6 


The zero position of the radii, for which Us and Uz are in phase, is 
defined by 


Z,=2(1- 


Ph 6° cS rLY 
Set eee ee 6 (63) 


The relations (58), (61), (62) and (63) determine the centers and radii 
of the circles of power with all the accuracy which has physical meaning 
for lines up to 300 miles, and without any necessity of using the hyper- 
bolic functions. 
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CONCLUSIONS 


We have shown that because of the errors in the constants of the 
line, generally there is no physical meaning in using more than two 
terms in the developments of Z, and Y,. Furthermore, because it is of 
no use to try to determine the magnitudes beyond a certain accuracy, 
in many cases it is perfectly correct to limit the expansions of Z, and 
Y, to the first terms, that is Z and Y/2, respectively. 

The definition of a line as a “short” (of: negligible capacitance) if 
the length is less than a given value, is not rational. It is necessary to 
specify the criterion employed. We have seen that a line is ‘‘short”’ if 


L < 20 miles by the condition of neglecting Js9 S 0.05 In, (V) 
L < 100 miles by the condition Uso/Uro > 0.98, (IV) 


L < 240 miles by the condition Iss/Uss5 exact 
within 10 per cent., (VI) 


for specified values of y, Ar/r, A(/w)/lw, AL/L. 
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NOTES FROM THE NATIONAL BUREAU OF STANDARDS.* 


DETECTION OF ISOTOPES BY TRACER MICROGRAPHY. 


A new method for the more effective tracing of radioactive isotopes 
in materials in which they have been intentionally introduced has been 
developed by L. Marton of the National Bureau of Standards with 
the cooperation of P. H. Abelson of the Department of Terrestrial Mag- 
netism, Carnegie Institution of Washington. In this procedure, by 
means of a magnetic focusing arrangement, the radiation given off by 
a radio-isotope within a sample material is made to form an image of 
the emitting surface upon a photographic plate. This image may then 
be used in studying the distribution and concentration of the radioactive 
element present in the sample. 

In many chemical, biological, biochemical, and other fields of re- 
search, there is growing application of the method of tracers, in which 
the isotope of a given element is used as an indicator to tag or label 
certain groups of atoms so that they may be distinguished from other 
atoms of the same kind. Identification of tracer elements is at present 
greatly facilitated through the use of radioactive isotopes, which, be- 
cause of recent developments in atomic energy, are now available in 
large quantities and are relatively easy to detect through their radia- 
tions. 

In the well-known method of radio autography a radio-isotope is 
introduced in a biological or other system, and the distribution of that 
particular element within the system is determined by bringing the 
sample in close contact with a photographic emulsion. This method 
lacks resolving power because, even in case of perfect contact of the 
sample with the emulsion, the circle of confusion from every point of 
emission is so great that details less than a tenth of a millimeter are very 
difficult or impossible to distinguish. 

In order to improve the resolution of this tracer method, it was 
decided to use electron optical image formation for determination of the 
distribution of a radioactive element within a given sample. This 
process, which may be called ‘‘tracer micrography,” is based on the 
emission of high speed electrons (beta rays) by many tracer elements 
and the use of magnetic lens elements to form an image on a suitable 
recording surface. 

In the absence of any means for correction of the chromatic aber- 
ration of electron optical lenses, the first micrographs were limited to 
those elements which emit electrons of uniform speed. After some at- 


* Communicated by the Director. 
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tempts with Columbium™, Yttrium®’, Strontium®, Strontium®’, and 
Protactinium**, Gallium*®’ was selected for the initial tests. Gallium 
chloride, prepared by chemical separation from zinc, was bombarded by 
heavy hydrogen nuclei in the cyclotron at the Carnegie Institution, and 
the solution was evaporated drop after drop on a }-inch tantalum disc. 
Radiation emitted from the surface of the disc, upon passing through a 
magnetic lens consisting of a small ironclad coil with Armco iron pole 
pieces, was brought to a focus upon a photographic film at a distance of 
about 33 inches. An image of the tantalum disc was thus obtained 
showing radioactive areas. The conditions were selected so that a 
linear magnification of 2 was obtained. 

For calibration of the instrument, the photographic film was re- 
placed by a Geiger counter, and the lens current necessary to produce a 
maximum number of counts in unit time was determined for radiations 
of varying velocities. This establishes the focusing current for a given 
type of radiation. 

_ In preliminary experiments with samples of different concentration 
and thickness of the radioactive layer, exposure times ranged from 2 to 
12 hours according to the age and concentration of the sample and the 
numerical aperture of the lens. It was found that micrographs with 
good definition were obtained consistently when the layer was suffici- 
ently thin to avoid considerable self-absorption. The best resolving 
power attained so far has been about 30 microns. 

The simplicity of this method, both in apparatus and technic, is one 
of its more important features. Vacuum requirements are very moder- 
ate, since the mean free path of the electrons is large in comparison with 
the apparatus dimensions, even at forepump pressure. 

Further improvements in tracer micrography are expected through 
after-dcceleration of the beta particles by means of an homogeneous 
electrostatic field. Such after-acceleration may well result in reduced 
exposure time and in better resolution due to a reduction in spherical 
aberration. A further reason for after-acceleration is that chromatic 
aberration, which is always present, even in sources emitting particles 
of uniform speed, can be markedly decreased if the accelerating potential 
is at least comparable in magnitude to the energy of the primary 
emission. 

Location of radioactive tracer elements is facilitated by a new elec- 
tron optical method, “tracer micrography.”’ 


TEMPERATURE DISTRIBUTION IN A TEST BUNGALOW. 


Uniformity of temperature throughout houses is a tacitly accepted 
American ideal of heating. It is probably never attained in practice, 
and departures from it depend upon the design and construction of the 
house and upon the characteristics of the heating system or device used. 
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Distribution of warmth in houses has usually been judged only qualita- 
tively by individual engineers on the basis of personal experience or 
observation. Quantitative data have been limited almost entirely to 
laboratory tests on parts of the system. In order to obtain data on the 
temperature distribution in an entire house, a series of heat distribution 
tests of several types of heating systems were conducted at the National 
Bureau of Standards by Richard S. Dill and Paul R. Achenbach. 

A full scale house in which complete heating systems can be installed 
was constructed at the Bureau. Designated as a test bungalow, it is 
similar in plan to house B described in F.H.A.’s Technical Bulletin No. 
4, “Principles of Planning Small House.”’ It has four rooms and bath, 
with a central hallway. The walls are conventional in construction 
consisting of 3-inch gypsum board on the inside and 2- by 4-inch 
studding with sheathing and lap siding on the outside, separated by a 
layer of building paper. All the windows are double hung with the ex- 
ception of that in the bathroom, and one of the two in the kitchen. 
The double floor of 1-inch pine includes building paper between the 
subfloor and finish floor. For the heat distribution tests the walls 
were not insulated but a 2-inch blanket of wood-fiber insulating material 
was laid over the ceiling. Weather-stripping is not provided for windows 
or doors. 

During the tests, data on temperature conditions inside the bungalow 
were recorded by various types of instruments. Heat-transfer coeffi- 
cients for floors, side walls and ceilings were measured by the Nicholls 
type of heat-flow meters fastened to the surfaces. Thermocouples 
located at various levels of the rooms, on top and bottom of the floor, 
on the inside and outside of the walls, and in the basement and attic 
air, gave a complete record of temperature conditions. Observations 
were made while heat was being supplied in turn by an experimental 
electric heater, an oil-burning warm-air furnace, a jacketed gas-firing 
space heater, a jacketed oil-fired space heater, a single gas-burning 
gravity floor furnace, two gas-burning gravity floor furnaces, a gas- 
burning floor furnace with forced circulation, an oil-burning gravity 
floor furnace, and a conventional gravity hot-water heating system. 

From the standpoint of comfort, the temperatures from the floor to 
5 feet above the floor are more significant than the temperatures at 
higher levels. The tests showed that the average temperature differ- 
ences produced in the test bungalow between the 2-inch level and the 
60-inch level, with continuous forced circulation of air through a plenum 
chamber, with a forced-circulation gas floor furnace, and with a gravity 
hot-water heating system, were less than 10°F. when the outside 
temperature was 32°F. The temperature differences produced in this 
same zone by all the other devices tested ranged from 14 to 18° F. at 
comparable outside temperatures. 

-The average horizontal temperature differences between rooms at 
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all levels of measurement ranged from 2° to 4° F. for the several floor 
furnaces and for the gravity hot-water system. For the electric warm. 
air furnace and the oil-burning furnace, when attached to a plenum 
chamber, these differences were nearly as small for the 2-, 30-, and 60- 
inch levels, but were considerably greater in the upper levels of the 
house. The temperature difference between rooms was greater for the 
space heaters in the living room than for the other heaters tested. This 
effect can probably be attributed to the fact that the rooms other than 
the one containing the heater are warmed by the overflow of heat from 
the living room through the doorways. Except for the space heaters, 
the average horizontal temperature differences from 2 inches to 6() 
inches above the floor did not increase with lower outside temperatures. 

Certain types of heaters may be suitable for houses of one construc- 
tion and not suitable for other constructions. For example, a high 
temperature near the ceiling would not cause much additional heat 
loss if the ceiling were well insulated, or an insulated floor might permit 
the use of heaters that do not deliver warmed air at or near the floor 
level. Furthermore, some types of heaters that provide comfort in 
areas where extremely low temperatures do not occur may not provide 
comfort in colder climates. 

Although standards of performance in the field of vertical or hori- 
zontal temperature differences in residences have not been established, 
the Bureau tests indicated that the air temperatures in the living zone 
of all rooms should be in the range from 65° to 80° F, for comfort when 
heated by the conventional types of heating systems. A floor tempera- 
ture of 60° F., when continued for 1 hour or more, was found to be too 
low for foot comfort; when a floor temperature of 65° F, was continued 
for the same period, discomfort was no longer apparent. It was further 
observed that air temperatures of approximately 85° F. around the 
level of the head became oppressive, especially if the air had an appreci- 
able velocity. 

The results obtained with the devices used in the test bungalow 
show need for further improvement in the design of house construction 
and heating systems in order to provide comfort in all rooms of base- 
mentless houses. A more uniform temperature can probably be ob- 
tained either by the use of more insulation in the house elements, by 
improving the heat-distribution systems, or by a combination of both. 

Entire heating systems were installed in this test bungalow, con- 
structed by the Bureau for studies of heat transfer phenomena, in order 
to determine the performance of various heating devices designed for 
use in low-cost housing. The bungalow is now enclosed with an in- 
sulated shell so that any desired temperature can be maintained outside 
of the house. 


Nore: The complete technical report on this project is contained in the National Bureau 
of Standards Building Materials and Structures Report BMS108, ‘‘Temperature Distribution 
in a Test Bungalow with Various Heating Devices.” 
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MERCURY-BALANCE METHOD FOR THE MEASUREMENT OF GEL STRENGTH. 


An improved method for more precise measurements of the tensile 
strength of gels’ has been developed by W. J. Hamer of the Bureau’s 
electrochemistry laboratory. This procedure, which gives the strength 
or “yield value” and resistance to deformation of the material under 
study, essentially involves the determination of the shearing force 
necessary to fracture the gel as given by the weight of mercury required 
to pull a standard disc from the gel in which it is embedded. 

In many applications of gelatin, glue, starch, agar, and various gums, 
it is important to know the tensile strength of the gel or jelly. This 
property was originally measured by comparison with a set of arbitrary 
standards on the basis of the resistance offered by the gel to finger pres- 
sure. Although this procedure was undoubtedly satisfactory for many 
commercial purposes, there was need for a more objective method. 
Such a method is very important, for example, in studies of gelatin 
coatings in photographic operations, particularly with the increasing 
use of flexible film supports instead of rigid glass bases. 

The new Bureau development represents objective method for gel- 
strength measurement. While this procedure was designed specifically 
for the study of the paste walls of dry cells, it may be applied to all types 
of gels and used in the study of food starches, in following small changes 
in starch pastes for use in adhesives, and in many other applications. 

The method involves the balancing of the cohesive forces of the gel 
against a gradually increasing quantity of mercury until the yield point 
is reached and the gel is broken. A brass disc is suspended in a beaker 
containing the gel by means of a wire attached to one beam of an analyt- 
ical balance. An adjustable platform over the balance pan supports 
the beaker. The gel is covered with liquid petrolatum to prevent “skin”’ 
formation. Mercury is added from a burette at a constant rate to a 
container on the other pan. Different flow rates are made possible by 
use of a removable capillary tip held to the burette by a rubber tube. 
The capillary tip is calibrated before use by weighing the mercury it 
delivers in a given time. As the mercury is added, the height of the 
platform is changed manually, so that the pointer of the balance is 
always kept at zero. A vertical scale indicating the height of the ad- 
justable platform is read at frequent regular intervals until the gel is 
broken. 

The weight of mercury required to produce a fracture in the gel 
divided by the area of the brass disc gives gel strength in grams per 
square centimeter. The rate of change in height of the adjustable 
platform gives an indication of the resistance the gel offers to deform- 
ation. A fast rate means that the gel is easily deformed, whereas a 

1 Details of this work are described in Bureau Research Paper RP1810, an improved 


method for measurement of gel strength and data on starch gels. W. J. Hamer, J. Research, 
NBS 39 (1947). 
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slow rate indicates that the gel resists deformation and is probably 
quite rigid. Although rigid gels usually offer high resistance to de- 
formation, this is not always true. 

This method has a precision greater than 0.5 per cent. Its limita. 
tions lie more in the reproducibility with which gels can be prepared 
than in errors inherent in the method. It may be applied to all types 
of gels, whether they are of the reversible (gelatin) or irreversible 
(starch) type. 

Gel-strength measurements are used in the control of several in- 
dustrial processes, such as the manufacture of adhesives, pectin jellies, 
starch sponges, and dextrins; in the preparation of foods; and to a more 
limited extent in the paper and textile industries. They are also fre- 
quently used in studies of the granular structure of starches and in 
determinations of the effects of acids, heat, aging, alcohols, polyhydrols, 
and other physical or chemical factors on the elasticity of starch and 
gelatin gels. 

The paste wall of most dry cells consists of a mixture of starch and 
flour in the form of an aqueous gel. Recently the Bureau has made 
studies of the paste-wall separator of dry cells of the Leclanche type. 
Gel-strength measurements were used to differentiate walls of different 
starches, to find out whether the paste reacted with the other constitu- 
ents of the cell, and to determine the effects of chemical and physical 
modifications on the strength of the paste wall. The results were cor- 
related with the shelf life of dry cells made with the various types of 
paste walls. 

The characteristics of the paste walls are often described in terms 
of gelling temperature of the paste or the variety and amount of starch 
and flour used in its preparation. Frequently the viscosity of the paste 
is also used for this purpose. For quantitative results, however, a 
more sensitive method is needed, preferably one that may be used with 
the gels themselves rather than the pastes from which they are formed. 
Numerous procedures have been proposed for objective measurements of 
gel strengths. Several of these are limited in application, because they 
depend on the breaking or compression of the surface of the gel and 
are subject to errors owing to skin formations of variable thickness. 
Furthermore, many methods give a measure of different characteristics 
of a gel depending upon the point of view of the experimenter. For 
example, they may indicate rigidity, elasticity, plastic viscosity, or 
resistance to some force such as cutting action or torsion. The new 
method developed at the Bureau gives consistent results, well adapted 
to use in the Bureau’s studies of the paste walls of dry cells. 
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THE FRANKLIN INSTITUTE. 


STATED MONTHLY MEETING, WEDNESDAY, DECEMBER 17, 1947. 


The stated monthly meeting of The Franklin Institute was held on Wednesday, December 
17, 1947, at 8:15 p.m. The President, Mr. Richard T. Nalle, presided over the meeting. The 
President announced that the minutes of the stated monthly meeting for October were printed 
in full in the November issue of the JouRNAL OF THE FRANKLIN INsTITUTE, and if there were 
no corrections or additions, the minutes would stand approved as read. The minutes were 
approved. 

The President then called on the Assistant Secretary, Dr. Frazer, for his report and 
announcements. Dr. Frazer announced that the following members had been elected during 
the month of November: 


Total membership as of November 30 


Dr. Frazer then made the announcement of our Christmas Week Lectures, which are given 
for young people through the James Mapes Dodge Foundation. It was announced that these 
lectures would be held on December 30th and 31st at 4:30 p.m., and that the speaker this year 
would be Mr. Vincent J. Schaefer, who would speak on “Scientific Fun in the Fields, the Moun- 
tains, and the Sky.” 

Dr. Frazer further announced that according to the By-Laws, Article IV, Section 5, 
nominations for the Board of Managers were to be made at this December meeting. The 
following members of the Board of Managers, to serve for a period of three years, were 
nominated: 

HiraM S. LUKENS 

Orvus J. MATTHEWS 
Morton GIBBONS-NEFF 
CHARLES S. REDDING 
GEORGE WHARTON PEPPER 
M. M. PRICE 

James S. ROGERS 
CLARENCE TOLAN, JR. 


The President then called for any nominations from the floor. Since no nominations were 
made, the President announced that the nominations for members of the Board of Managers 
were Closed. 

The President then presented the speaker of the evening, Dr. Donald H. Andrews, Pro- 
fessor of Chemistry at the Johns Hopkins University, in Baltimore, Maryland. Dr. Andrews 
spoke on “Radio and Infra-red Detection by Superconductivity.”” Dr. Andrews demonstrated 
a new type of bolometer, which was developed during the recent war at the Johns Hopkins 
University. The bolometer was used for detecting infra-red rays. It consists of a small 
ribbon of columbium nitride cemented to the top of a copper block 1 cm. square and 1 cm. high. 
When cooled to —258.5° C., the ribbon becomes partially superconducting and detects both 
infra-red and radio waves with a sensitivity in many ways superior to any other known device. 
Recordings of various kinds of signals received by the bolometer were demonstrated both with 
an oscilloscope and a loud speaker. 

After Dr. Andrews’ very edifying lecture, the President adjourned the meeting with a rising 
vote of thanks to the speaker. 
JOHN FRAZER, Assistant Secretary. 
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New MEMBERS. 


NEW MEMBERS OF THE FRANKLIN INSTITUTE 
ELECTED JUNE 1, 1947 TO DECEMBER 1, 1947. 


Mr. Sherman Mills Fairchild 


SUSTAINING MEMBERS. 


Mr. Eugene J. Houdry 


ACTIVE MEMBERS WITH FAMILY PRIVILEGES. 


Mr. Pierre Aureli 

Mr. N. Newlin Baily 

Mr. Robert W. Barker 
Mr. Mayer I. Blum 

Mr. Elmer J. Boetefuer 
Mr. Edward G. Budd, Jr. 
Mr. Emil Cohn, Jr. 

Mr. John M. Diem 

Mr. Moses Ehrlich 

Mr. Joseph H. Gillies 


Mr. William Loughridge 
Aiken 

Mr. Edward J. Bailey 

Major General Milton G. 
Baker 

Mr. Henry S. Bamford 

Mr. Naaman F. Barr 


Mr. William L. Baumner, Jr. 


Mr. Charles W. Berg 
Mr. Sanford Bliwise 
Mr. Lou Block 

Mr. M. A, Bordman 
Mr. Harry S. Boud 
Mr. John E. Bower 


Mr. Wiliiam H. Brearley, Jr. 


Mr. W. Ronald Briggs 
Mr. Harold W. Brown 
Mr. James M. Castle 
Mr. George A. Clickenger 
Mr. Frank G. Conrad, Jr. 
Mr. Edward F. Cooper 
Mr. Edward P. Cornely 
Mr. George J. Crits 

Mr. Joseph Croskey 

Mr. E. B. Curdts 

Mr. Edwin A. Dages 

Mr. Joseph De Frenes 
Mr. William J. Devinney 
Mr. Felix Di Arenzo 

Mr. Ralph Earle 

Mr. A. J. Ellis 

Mr. Jay Emanuel 

Paul J. Ernst, Ph.D. 


Mr. Charles H. Godschall 
Mr. Harry C. Hambridge 
Mr. Harry R.!Hirshorn 
Mr. David H. Hopkins 
Mr. Laurence T. Howell 
Mr. Edward Kohlhepp 
Mr. George Lipsky 

Mr. John S. Malick 

Mr. John Cochrane Martin 
Mr. John F. McNelis 

G. Herbert Moffses, M.D. 


ACTIVE MEMBERS. 


Mr. Woodfin R. Fagge 

Mr. William F. Famille 

Mr. Mort F. Farr 

Mr. Myer Feinstein 

Mrs. Jeanne Ferber 

Mr. Alex C. Fergusson 

Mr. Harry B. Fertik 

Thomas Fitz-Hugh, Jr., 
M.D. 

Mr. Lawrence S. Ford 

Mr. John J. Fox 

Mr. Benjamin H. Freeman 

Mr. Charles H. Fryburg 

Mr. Edward E. Gallob 

Mr. Arthur S. Gleason 

Mr. Richard A. Griefen 

Mr. H. K. Groff 

Mr. William P. Groves, Jr. 

Mr. Amil S, Gumula 

Mr. Fred C. Haab 

Mr. Joseph J. Hagan 

Mr. Wilson C. Hanline 

Mr. Paul V. Hannas 

Mr. Abdullah Feyzi Hamdi 

Mr. W. J. Heggie 

Mr. Robert F. Herr 

Gerald E. Herrnstadt, Ph.D. 

Mr. Emil Hoffman 

Mr. William P. Hood 

Mr. William Hopkin 

Mr. Wayne L. Hopkins 

Mr. John J. Howard 

Mr. A. D. Howry 


Mr. Harry E. Mousley 
Mr. Robert’ Nagel 

Mr. Gordon Null 

Mr. Marvin H. Pond 

Mr. Howard C. Pruyn 
Mr. Harry Rubin 

Mr. Theodore A. Seraphin 
Mr. Samuel A. Weston 

O. Wilson Winters, D.D.S, 
Mr. Carl A. Wolf 


Mr. Robert V. Hudson 

Mr. Joseph Gray Jackson 

Mr. Bernard Jacobs 

Mr. C. Brooke Jones 

Mr. Morris Katz 

Mr. John Keefe, Jr. 

Mr. John Kerbeck 

Mr. Ellis Knowles 

Mr. William Kolb 

Mr. Alfred Theodore Korn- 
field 

Mr. Ellis M. Landis 

Mr. Henri Lauer 

Miss Dorothy K. Lawrence 

Mr. Hurley Lazarus 

Mr. Alfred Leith 

Mr. Adam Leva 

Norton M. Levin, D.O. 

Mr. Willard F. Levy 

Mr. Joseph F. Lewis 

Mr. John Franklin Long 

Mr. Thomas C. Longstreth 

Mr. James Luckman 

Mr. Lionel P. Maggenti 

Mr. Charles W. Maloney 

Mr. John Mallo 

Mr. Thomas I. Manly 

Mr. Samuel A. Mann 

Mr. Vernon L. Marquart 

Mr. Emil A. Mathias 

Mr. T. E. McBride _ 

Mr. Laurence T. McCurdy 

Mr. Sterling R. Mensch 


Harold | 
Mr. Mo 


Mr. Gec 
Mr. Car 
Mr. Joh 


Mr. Cha 
Miss Eli 
Mr. C. F 
Mr. Geo 
"44 
Mr. Juli 
Mr. Fre 
41 
Mr. Rot 
Mr. Mor 
Mr. Wal 
Mr. Cha 
Mr. Will 
Mr. Joh 
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Mr. C. B. Miller 

Mr. Raymond Millman 

Mr. Harry W. Morgan 

Mr. Robert V. Morris 

Mr. Richard J. Mott 

Mr. Arthur Niessen 

Mr. Charles H. Nonamaker 
Mr. M. John O'Donoghue 
Mr. Gordon Palmer, Jr. 
Mr. Daniel C. Parr 

Mr. Edward Paul 

Mr. Richard D. Pomerantz 
Mr. Clarence H. Porter 
Mr. F. B. Powers 

Mr. Harold C. Rahn 

Mr. Leonard Rodenhausen 
Mr. Frank B. Russell 

Mr. Herman Schaevitz 
Harold G. Scheie, M.D. 
Mr. Morris Schiff 


Mr. George W. Alcock 
Mr. Carleton Craig 
Mr. John R. Heizmann 


Mr. Charles S. Beach '01 

Miss Elizabeth C. Biddle '09 

Mr. C. H. Bierbaum '97 

Mr. George Imlay Bodine, Jr. 
"44 

Mr. Julius Boekel '84 

Mr. Franklin S. Chambers 
"41 

Mr. Robert Cherry, Jr. '37 

Mr. Morris L. Clothier '26 

Mr. Walter T. Dwyer '46 

Mr. Charles J. Eisenlohr '36 

Mr. William M. Elkins ’42 

Mr. John J. Foulkrod, Jr. ’36 


Mr. Frederick J. Schneider 
Mr. Howard I. Scherr 

Mr Gerry W. Schulz 
MrJ4Richard Schweder 
Mr.4Robert D. Scott 

Mr. Louis Segal 

Mr. John B. Shields 

Mr. Chao Tze Shih 

Dr. Alfred R. Smith 

Mr. Cortland Gray Smith 
Mr. Edward W. Smith, Jr. 
Mr. Perry C. Smith 

Mr. Albert Soffa 

Mr. Robert Solts, Jr. 

Mr. Ulrich Sontheimer 
Mr. James C. Stambaugh 
Mr. Albert D. Stave 

Mr. Saul H. Steel, Jr. 

Dr. John A. Stevenson 
Mr. Charles M. G. Stewart 


ACTIVE NON-RESIDENT MEMBERS. 


Mr. Arthur T. Hinckley 
Mr. Mortimer Francis Nahon 
Mr. Guy F. Rolland 


NECROLOGY. 


Mr. William H. Frey '43 

Mr. Anthony H. Geuting 36 

Mr. John Elgin Gilmer 43 

Mr. Joseph Handler '43 

Mr. Frederick A. Healy,’36 

Boyden W. Kowalski, M.D. 
44 

Mr. John M. Krupp ’43 

Mr. James Barton Longacre 
"38 

Mr. Sylvester A. Mahan '46 

Mr. James T. J. Mellon '42 

Samuel Nicholas, M.D. '44 

Max Planck, Ph.D. '27 

Mr. J. Benton Porter '41 


COMMITTEE ON SCIENCE AND THE ARTS. 


Mr. E, Lee Strock 

Mr. Lisle J. Sykes 

Mr. John J. Townsend 
Mr. John W. Temple 
Mr. John M. Tinker 
Mr. Irving S. Towsley 
Mr. Wilmer S. Trinkle 
Mr. William A. Tucker 
Mr. Russel R. Tull 

Mr. Edward S. Turner 
Mr. Ralph Vezin 

Mr. Alfred M. Whitney 
Mr. Edward J. Williams 
Mr. R. M. Wilson 

Mr. Samuel S. Wolfman 
Mr. R. &. Worden 
Carroll S. Wright, M.D. 
Mr. Joseph A. Wurster, Jr. 
Mr. Samuel M. Zollers 


Mr. P. A. Shields 
Mr. George H. Woodard 


Mr. Paul C. Richter '46 

Mr. Isaa¢ B. Ritter '42 

Mr. Charles E. Schmidt ’43 

J. E. Shrader, Ph.D. '22 

Mr. J. Stogdell Stokes ’36 

Mr. S. M. Swaab '46 

Mr. Henry F. Wageman °46 

Miss Hepsey Norris Wells '34 

Mr. William H. B. Whitall 
36 

Thomas L. Wilcox, D.D.S., 
"45 

Major Richard R. Wright, 
Sr. '46 


COMMITTEE ON SCIENCE AND THE ARTS. 


(Abstract of #roceedings of Stated Meeting held Wednesday, November 12, 1947.) 


HALL OF THE COMMITTEE, 
PHILADELPHIA, NOVEMBER 12, 1947, 


Mr. WALTER C. WAGNER in the Chair. 
The following report was presented for final action: 
No. 3186: Henderson Medal. 


This report recommended the award of a George R. Henderson Medal to Charles Dun- 
canson Young, of Whitford, Pennsylvania, ‘In consideration of his contribution to the 
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scientific advancement of the steam locomotive which has resulted in improving the reliability 
and efficiency and reducing the cost of steam locomotives thereby producing a more effective 
transportation unit.” 


(Abstract of Proceedings of Stated Meeting held Wednesday, December 10, 1947.) 


HALL OF THE COMMITTEE, 
PHILADELPHIA, DECEMBER 10, 1947. 
Mr. WALTER C, WAGNER in the Chair. 
The following report was presented for final action: 
No. 3185: Electrical Gun Director. 

This report recommended the award of a Howard N. Potts Medal each to David Bigelow 
Parkinson, of Maplewood, New Jersey, and Clarence Anding Lovell, of Summit, New Jersey, 
“In consideration of their combined contributions, both to the theoretical and practical design 
of the Electrical Gun Director, one of the outstanding pieces of equipment developed during 
the war period, which development has further contributed largely to the theory and practical 
application of servo mechanisms and smoothing filters in general.” 

JoHN FRAZER, 
Secretary to Commmittee. 


LIBRARY NOTES. 


The Committee on Library desires to add to the collections any technical works that 
members would wish to contribute. Contributions will be gratefully acknowledged and placed 
in the library. Duplicates received will be transferred to other libraries as gifts of the donor. 

Photostat Service. Photostat prints of any material in the collections can be supplied 
on request. Orders received in the morning are filled the same day. The average cost for 
a print 9 X 14 inches is thirty-five cents. 

The Library and reading room are open on Mondays, Tuesdays, Fridays and Saturdays 
from 9 A.M. until 5 P.m., Wednesdays and Thursdays from 2 p.m. until 10 P.M. 


RECENT ADDITIONS. 
BACTERIOLOGY. 
HinsHELWooD, C. N. The Chemical Kinetics of the Bacterial Cell. 1946. 
CHEMISTRY AND CHEMICAL TECHNOLOGY. 


CHEMICAL RuBBER CoMPANY. Handbook of Chemistry and Physics. Thirtieth Edition. 

CoLiorp Science. A Symposium. 1947, 

Dor&e, CHARLES. Methods of Cellulose Chemistry Second Edition. 1947. 

EAsTMAN, E. D., AnD G. K. RoLLEFsSON. Physical Chemistry. 1947. 

EITEL, WILHELM. Physikalische Chemie der Silikate. 1941. 

HEATON, Noe. Outlines of Paint Technology. Third Edition. 1947. 

KnaccGs, NELtson S. Adventures in Man’s First Plastics. 1947. 

Ko.LtHorr, IsAAc, AND V. STENGER. Volumetric Analysis. Second Edition. Volume 2. 
1947. 

MANTELL, CHARLES L. The Water-soluble Gums. 1947. 

Preters, H. A., AND D. W. VAN KREVELEN. The wet Purification of Coal Gas and Similar 
Gases. 1946. 

PricgE, CHARLES C. Mechanics of Reactions at Carbon-Carbon Double Bond. 1946. 


REMINGTON, JOHN STEWART. Pigments. 1944. 
ScumipT, J.,AND H. G. Rute. A textbook of Organic Chemistry. Fifth Edition. 1947. 


SoMMER, Huco Henry. The Theory and Practice of Ice Cream Making. Fifth Edition. 
1946. 
THIENE, HERMANN. Glas. Volumes 1 and 2. 1931 and 1939. 
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ELECTRICITY AND ELECTRICAL ENGINEERING. 


Barrows, E. J. Industrial High Voltage. Distribution and Public Supply. 1946. 
CuLVER, CHARLES AARON. Theory and Applications of Electricity and Magnetism. 1947. 
L’HERMITE, RoBERT. L’experience et les Théories Nouvelles en Resistance des Matériaux. 


1945. 
Laurson, P. G., A. MERRILL, AND W. JAMEs. Elements of Mechanism. Second Edition. 


1947. 
ORCHARD, FREDERICK C. Electricity in the Building Industry. 1946. 
SLATER, J. C., AND N. H. FRANK. Mechanics. 1947. 
WaLL, T. F. ELectricat ENGINEERING. 1947. 


INDUSTRIAL MANAGEMENT. 
TREWMAN, Harry F. Mechanical Inspection. 1946. 


MATHEMATICS. 
Hoet, Paut G. Introduction to Mathematical Statistics. 1947. 
MARSHALL, W. R., AND R. L. Picrorp. Application of Differential Equations to Chemical 
Engineering Problems. 1947. 
METALLURGY. 
BREARSLEY, HARRY. Talks About Steelmaking. 1946. 
CuesTeR, J. H., AnD M. W. Turinc. Influence of Port Design on Open-Hearth Furnace 


Flames. 1946. 
Jenkins, Ivor. Controlled Atmospheres for the Heat Treatment of Metals. 1946. 
Netson, J. H. anD H. Smmman. The Application of Radiant Heat to Metal Finishing. 1945. 


Woop, V. N. Metallurgical Materials. 1946. 
OPTICS. 
ToLANSKY,S. High Resolution Spectroscopy. 1947. 
PHOTOGRAPHY. 


DisERENS, CHARLES. Handbuch der Photographie. Volume 1. No date. 
DISERENS, CHARLES. ‘Traite de Photographie. Volume 1. 1943. 
Morcan, W. D., H. M. Lester AND OTHERS. Graphic Graflex Photography Eighth Edition. 
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PHYSICS. 


MICHENER, WILLIAM H. Physics for Students of Science and Engineering. 1947. 


Morsg, F. T. Elements of Applied Energy. 1947 
SLATER, J. C., AND N. H. FRANK. Mechanics. 1947. 
Weser, R. L., M. W. Waite, AND K. V. MANNING. College Technical Physics. 1947 


RADIO. 
Harrison, ARTHUR EpwWARD. Klystron Tubes. 1947. 
SCIENCE. 


HALDONE, J. B.S. Science Advances. 1947. 
Heatu, H. F., AND A. L. HEATHERINGTON. [Industrial Research and Development in the 
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NOTES FROM THE BIOCHEMICAL RESEARCH FOUNDATION. 


Enzymatic Hydrolysis of Desoxyribonucleic Acid Prepared with the 
Use of Strong Sodium Hydroxide.—CuarLes A. ZITTLE. It was 
previously reported (1) that desoxyribonucleic acid (DNA) prepared 
under mild conditions (method of Hammarsten) (2) was readily hydrol- 
yzed by specific nuclease so that it became soluble in HCl, whereas 
DNA prepared with the use of hot 1.25 N NaOH (method of Levene) 
(3) was not hydrolyzed by the nuclease. Further study has shown 
that under the proper conditions the alkali-treated DNA can be slowly 
hydrolyzed by the nuclease at a rate about 20 per cent. of that of DNA 
(Hammarsten). The previously reported erroneus conclusion that 
DNA (Levene) could not be hydrolyzed was in consequence of this 
low rate of hydrolysis and the instability of the enzyme under the con- 
ditions used. 

Experimental. 


Preparation of the enzymes and substrates and the manometric 
and other procedures have been described (4). In initial experiments 
with DNA (Hammarsten) this nucleic acid was contemplated as a sub- 
strate for the phosphoesterase from calf intestinal mucosa (5), and for 
this purpose it was treated with specific nuclease since hydrolysis of 
DNA by the phosphoesterase is only rapid and complete after it has 
been acted on by specific nuclease (4). For preparing the nuclease- 
treated DNA the nuclease was obtained in the form of a dried am- 
monium sulfate precipitate and weighed amounts of this material were 
added directly to the nucleic acid soultions. The above observations 
about hydrolysis were made under these conditions. In subsequent 
attempts to use dilute solutions of the nuclease in manometric experi- 
ments it was found to be very unstable, but could be stabilized by the 
addition of neopeptone following observations of McCarty (6). 

A stable nuclease solution was prepared as follows: 300 mg. were 
dissolved in 20 cc. of 0.0005 N H,SO, and dialyzed against the same to 
remove the ammonium sulfate; 6.0 cc. of this solution were mixed with 
4.0 cc. of neutralized neopeptone, 1.3 cc. of 0.5 M NaHCO, and 0.7 cc. 
of 0.3 M magnesium chloride. This solution, prepared from the dia- 
lyzed enzyme for each experiment, was used in the experiments to be 
described. 

The initial rates of hydrolysis of the two nucleic acids were deter- 
mined manometrically. In each experiment 2.0 cc. of 0.6 per cent. 
nucleic acid together with 0.2 cc. of 0.5 M NaHCOs, variable amounts 
of enzyme in the side arm of Warburg flasks, and water to make the 
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final volume 3.5 cc. were used. In these experiments the initial rates 
of hydrolysis were: for DNA (Hammarsten), 91.2 cmm./10 min./1.0 cc. 
of nuclease; for DNA (Levene), 17.0. 

Precipitates obtained by acidifying solutions of DNA (Levene) 
with HCl turn a light brick red on standing for some time; initially 
the precipitates are white. This observation and the slow rate of 
hydrolysis of this nucleic acid suggested that an inhibitory substance 
might be present from the reagents (colloidal iron, for example) used 
in its preparation. Hydrolysis of the two preparations of nucleic acid 
in the same flask (followed manometrically: 6.0 mg. of each nucleic 
acid, 0.5 cc. of enzyme), however, was equal to the sum of the two taken 
alone, indicating that the alkali-treated nucleic acid preparation was not 
inhibitory to hydrolysis of the other. With the same thought in mind 
the alkali-treated nucleic acid was dialyzed for 24 hours, precipitated 
with HCI (final concentration, 0.24 N; pH 1.2) and dried with acetone. 
Fifty-eight per cent. of the starting material were recovered. The 
enzymatic hydrolysis of this material (followed manometrically: 3.0 to 
9.0 mg. of nucleic acid; 0.5 cc. enzyme), however, was the same as that 
of the original. 

For complete hydrolysis of the alkali-treated nucleic acid at 23°, 
measured by solubility in 0.25 N HCl, the following procedure was 
adequate. Nine parts of 0.6 per cent. solution of the nucleic acid were 
treated with successive portions of the enzyme solution: with two parts 
initially, with two parts at the end of 4 hours, and with two parts at the 
end of 7 hours. The last portion of enzyme was permitted to act for 18 
hours with the solution saturated with chloroform. Although precipi- 
tation with HCI was now negative, precipitation with the uranium re- 
agent had remained unchanged. 

Phosphoesterase was permitted, as was done previously (1), to act 
on the above digest without and with silver nitrate present, the silver 
to inhibit the adenosine deaminase in these preparations. With this 
complete digest the results are comparable to those obtained with digest 
of DNA (Hammarsten); with silver present more CO, is released due 
to the inhibition of the deaminase and consequent production of am- 
monia. This is in contrast to previous results with DNA (Levene). 
From the previous results it was difficult to decide whether the DNA 
had been deaminated by the strong alkali used in its preparation or that 
the phosphoesterase had not liberated adenosine. From the present 
data it appears that the latter explanation was the correct one. 


Discussion. 


The previous observation (1) that DNA (Levene) was not hydro- 
lyzed was made because of the instability of the enzyme under the con- 
ditions used and also because of the low rate of hydrolysis of this prepa- 
tion of nucleic acid. The DNA (Hammarsten) was rapidly hydrolyzed 
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and hydrolysis was complete before the activity of the enzyme was 
lost, but, the DNA (Levene) showed little or no hydrolysis during the 
same period of time. The results of the present study indicate that 
adenosine deaminase can act on both of these hydrolysates and ther 
is no reason for believing that the adenine radical is involved in the low 
rate of hydrolysis. 

Summary. 


In contradiction to results previously reported, it has been found 
that DNA prepared by a method which utilizes strong NaOH is acted 
on by the specific nuclease, although at a rate only 20 per cent. of that 
of DNA prepared by a milder method. The nature of the changes 
responsible for the low rate of hydrolysis of the alkali-treated DNA 
are not known. 
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BOOK REVIEWS. 


Arr CONDITIONING AND ELEMENTS OF REFRIGERATION, by Samuel P. Brown. First Edition. 
644 pages, drawings and illustrations, 15 X 23 cms. New York and London, McGraw- 
Hill Book Co., Inc. 1947. Price $6.00. 


“The objective of this book is to provide a complete, practical, working textbook and hand- 
book suitable for the study by any interested person and to those working in design, installation, 
operaticn of heating, ventilating, and other mechanical equipment of buildings.” The selection 
of equipment, automatic controls, and design of air ducts and liquid piping systems are covered, 
as well as refrigeration theory. 

The tables necessary for elementary design work are included, although complex theory 
and higher mathematics have been eliminated. 

Throughout the text, illustrative practical problems are presented, which are then answered 
in the appendix. There are over 230 graphs and illustrations throughout the book. 

Henry E. MICHAEL. 


FUNDAMENTALS OF ELECTRICITY AND MAGNETISM, by Leonard B. Loeb. Third Edition. 
669 pages, drawings, 15 X 24 cms. New York, John Wiley & Sons, Inc., 1947. Price 
$6.00. 


This textbook, originally written to fill a specific need at the University of ‘California, 
presents the subject, in its third edition, as a consistent and related whole, rather than isolated 
and dogmatically stated topics. 

The historical survey of previous editions is replaced by a brief historical outline and the 
space thus gained is devoted to the enlargement of the section on atomic structure. 

The book begins with elementary magnetism and direct current electricity followed by 
the study of static electricity. Thec.g.s. is utilized, since it is felt by the author that electricity, 
basically so closely associated with electronic and atomic phenomena, profits more from its use 
rather than the m.k.s. system, which has its advantages in application to gross mechanics or 
engineering physics. 

Henry E, MICHAEL. 


TABLES OF INTEGRALS AND OTHER MATHEMATICAL Data, by Herbert Bristol Dwight. Re- 
vised edition. 250 pages, tables, 14 X 21 cms. New York, The Macmillan Co., 1947. 
Price $2.50. 


This book contains a variety of mathematical data in formula and tables including 
algebraic, trigonometric, inverse trigonometric, exponential, logarithmic, hyperbolic, inverse 
hyperbolic, elliptic, and Bessel functions. There are also probability integrals, surface zonal 
harmonics, definite integrals and differential equations. An appendix gives tables of numerical 
values and a list of 62 references referred to in the text where the derivation of the results is 
siven and where further results may be found. At the end of the book there is also a subject 
index giving more details on material location so that for instance gamma function or Bernoulli’s 
tumbers may be found quickly. This second addition of this work is a convenient tool for 


those who use mathematics to any extent. 
R. H. OpPERMANN. 


Tue CHEMISTRY OF COMMERCIAL Puastics, by Reginald L. Wakeman. 836 pages, tables and 
illustrations, 15 X 24 cms. New York, Reinhold Publishing Corporation, 1947. Price 
$10.00. 


Twenty years ago it was not uncommon to hear it said that all pioneering in organic chem- 
stry was over—organic chemistry was no longer a virile science for no discoveries of funda- 
mental importance had been made by its devotes since the brilliant éclat of the aniline dye 
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industry. But the foundations were being laid for a new glorious period for organic chemistry 
the Age of Plastics. Plastics, as we know them today, cover a large group of organic sub. 
stances, either wholly or in part synthetic, which can be molded into coherent solid articles tha; 
retain their shape indefinitely at room temperature. This book is an analysis of the fund oj 
information on the subject which correlates the scientific data with actual present industria| 
practice. 

After an introductory chapter and a reference to the history of the plastics industry the 
work examines in a general way some of the relationships which exist between the various ray 
materials and the plastics made from them. This leads to treatments on general chemical 
principles of resinification and ways in which synthetic resins and plastics are manipulated jn 
making diverse products. Thereafter consideration is given separately to specific subjects 
such as phenolic, urea and melamine plastics, aniline-formaldehyde resins, polyamides-nylon, 
coumarone-indene resins, polyvinylidene chloride, polystyrene, rubber derivatives, natural 
and synthetic elastomers, cellulose, lignin and protein plastics, polymeric organosilicon-oxygen 
compounds. These represent the most important plastics of commerce. There have been 
developed many other synthetics of more or less resinous character and are marketed for special 
purposes in relatively small volumes. Some of these are taken up in the last chapter together 
with a few which have been studied extensively in industrial laboratories but which still have 
not appeared on the market. In the back of the book there is an author and a large well ar. 
ranged subject index, 

This book is the result of an exhaustive survey of the literature and an intimate knowledge 
of the industry. The many illustrations, tables, curves, diagrams coupled with skill in presen- 
tation and a progression to present practice pointing to the future make this book a valued 
contribution to the literature as a text and reference work. R. H. Opperman. 


TEXTILE Fipers by J. M. Matthews and H. R. Mauersberger. Fifth Edition. 1133 pages, 
tables and illustrations, 15 X 24 cms. New York, John Wiley & Sons, Inc., 1947. Price 
$12.50. 

The history of fiber development has been a strange pattern of trial and error ever since 
fibers were first used for making apparel. Climate had much to do with it since wool, hair 
and fur fiber had great and continual development in cold climates and in hot climates develop- 
ment was first found in vegetable fibers. With the invention of machinery, fiber development 
was more rapid. Some fibers were delayed in development because of man’s inability to dye 
them satisfactorily and to degum with sufficient thoroughness. Then came the need, partic- 
ularly of the totalitarian nations, of producing ersatz materials to take the place of natural 
fibers. New methods, new products, new machines, and new techniques continually displace 
the old. This volume on fibers depicts the progress that has been made in textile fibers, their 
economic utilization, and technical development. 

Now in its fifth edition, the book has 1095 pages of text, written for the student in colleges 
and technical school and as a reference for those who deal with the sale, production, processing, 
and conversion into fabric of all the different fibers now used in the textile and clothing in- 
dustries. A suitable introduction leads into a general detail of the essential physical, micro- 
scopic, and chemical properties which make fiber a valuable textile material. Later, these 
properties are given considerable space with regard to cotton, the bast fibers, structural or 
hard vegetable fibers, wool, specialty hair fibers, silk, rayon, new synthetic textile fibers, and 
material or inorganic fibers and filaments. A complete coverage is given including the eco- 
nomics, production, cultivation, statistics, prices, dyeing, spinning, etc. A considerable 
section is devoted to fiber identification methods, the quantitative determination of fibers 
present in a product, and fiber testing methods. There is a bibliography at the end of eacl 
of the 24 chapters and in the back of the book there is a large comprehensive subject index. 

This monumental work is the sum of the work of forty-six contributors, collaborators 
and reviewers whose work was supervised in selection, revision and arrangement by an editor 
with much experience in the industry. It is of practical value today and points to fields 0! 
research and promise for the future. R. H. Opraemann 
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SERVOMECHANISM FUNDAMENTALS, by Henri Laurer, Robert Lesnick, and Leslie E. Matson. 
First Edition. 277 pages, drawings, 16 X 23 cms. New York and London, McGraw- 
Hill Book Co., Inc., 1947. Price $3.50. 

The term“servomechanism”’ is generally used to denote a particular class of automatic 
control system where the controlling signals of the system change fast and frequently. This 
type of system is different from the “regulator’’ in that in the case of the regulator the con- 
trolling signals do not change or change slowly and the primary problem is to prevent ex- 
traneous causes from affecting the function controlled by the system. Servomechanisms are 
used to steer ships, control airplanes, automatically tune radio receivers, position guns, etc. 
Regulators control temperatures, maintain liquid levels etc. This newer servomechanism 
principle has the outstanding features that the controlled and controlling elements of the sys- 
tem may be widely separated, that the work performed may be many times greater than re- 
quired for control operation and the often complicated processes can be made to take place 
according to predetermined standards. 

This book is designed for practicing engineers as an introduction to the principles under- 
lying the theory of servomechanisms. As such it starts with the elements and proceeds to the 
more complex with sympathy to the needs of those who require refreshing or review coverage 
on many forgotten subjdcts. In detail, the work begins with defining the concepts of control 
of open and closed control systems and the means employed for accomplishing control opera- 
tions. Simple practical examples are described, enough to give an idea of duties performed 
which could hardly be accomplished by other means. The servo system follow up link; a 
differential device which produces an indication of the magnitude and direction of the error 
between the input and output members of the system, is next taken up. 

Because servo-control systems involve a knowledge of both mechanical and electrical 
engineering principles, the authors have inserted at this point a chapter on fundamentals of 
mechanics and electricity as a review covering only the most elementary concepts. Then pro- 
ceeding onward the'book treats on principles of the transient analysis method of elementary 
servomechanisms. This is divided as follows: analysis of servomechanisms with viscous out- 
put damping, with error rate damping, with combined viscous output damping and error-rate 
damping, with integral control, and an enlargement in greater detail devoted to practical 
means for obtaining error rate damping, and a discussion of error rate stabilization networks. 
There is an introduction to the transfer function analysis of servomechanisms which is a fre- 
quency analysis and the relation of this to the transient analysis is indicated. The last chapter 
of the book is devoted to typical design calculations and general consideration. At the end 
there is a subject index. 

The practicing engineer who has a need of familiarity with servomechanisms will find in 


this book a useful tool. 
R. H,. OpPERMANN. 


PETROLEUM PropUCTION, by Park J. Jones. Volume III. 271 pages, 16 X 24 cms., tables 
and drawings. New York, Reinhold Publishing Corporation, 1947. Price $5.00. 

This is Volume III of a series of books and is limited to oil production by water. This 
excludes primary gas caps and oil production by gas. The book is divided into two parts, the 
first dealing with foundations. Essentially in the form of a text, formulas, curves, and charts 
are given with discussions on reservoir space in barrels, the fluid factor and its effects, oil pro- 
duction by expansion, well producing capacity and the factors effecting it, injection of water, 
the optima for reservoirs involving economical as well as physical factors, and the migration of 
oil. Part two of the book is on applications and its coverage includes discussions on radial 
reservoirs, two pay intervals, impermeable wedges, saturated oil, linear reservoirs, elongated 
reservoirs, and bottom water. A subject index is in the back. 

The presentation shows a practical application of facts and figures relating to this type of 


petroleum production. 
R. H. OPPERMANN. 


PUBLICATIONS RECEIVED. 


The Chemistry and Technology of Waxes, by Albin H. Warth. 519 pages, drawings, table 
and illustrations, 15 X 23 cms. New York, Reinhold Publishing Corporation, 1947. Price 
$10.00. 

Architectural Construction, by Theodore Crane. 414 pages, drawings and illustrations 
16 X 25cms. New York, John Wiley & Sons, Inc., 1947. Price $6.00. 

Field Practice, by Elwyn E. Seelye. A Data Book'for Civil Engineers. 306 pages, tables 
and drawings and illustrations, 13 X 21 cms. New York, John Wiley & Sons, Inc., 1947 
Price $4.50. 

Saga in Steel and Concrete, by Kenneth Bjork. 504 pages, 16 X 24 cms. Northfield, 
Norwegian-American Historical Association, 1947. Price $4.00. 

Industrial Health Engineering, by Allen D. Brandt. 395 pages, drawings, tables and illus. 
trations, 15 X 24cms. New York, John Wiley & Sons, Inc., 1947. Price $6.00. 

Flight Engineering and Cruise Control, by Harris G. Moe. 209 pages, drawings and illus- 
trations, 15 X 23 cms. New York, John Wiley & Sons, Inc., 1947. Price $4.00. 

Very High-Frequency Techniques, compiled by the staff of the Radio Research Laboratory, 
Harvard University. VolumesI and II. First Edition. 554 pages, and 1057 pages, drawings 
and illustrations, 16 X 24 cms. New York, McGraw-Hill Book Co., 1947. Price $14.00, 
two volumes. 

The American Annual of Photography 1948, edited by Frank R. Fraprie and Franklin | 
Jordan. 216 pages, illustrations, 18 X 25 cms. Boston, American Photographic Publishing 
Co., 1947. Price $2.00. 


Chemical Engineering Fundamentals, by Chalmer G. Kirkkride. 419 pages, tables and 
drawings, 15 X 24 cms. New York, McGraw-Hill Book Company, Inc., 1947. Price $5.00. 
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CURRENT TOPICS. 


Detecting Impurities in Uranium. (Electrical Engineering, Vol. 66, No. 
7.) Thirty-three impurity elements, some in concentrations as low as one 
part in 20 million, have been detected and estimated through a modified spec- 
trographic method developed at the National Bureau of Standards and used 
since 1942 in the Manhattan Project for analysis of uranium and its compounds. 

Impurities differ in their effects on the nuclear chain reaction. Some ele- 
ments such as boron and cadmium may interfere if present in concentrations as 
low as a few tenths of a part per million, and many other elements should not 
exceed a few parts per million. Rapid, sensitive, and accurate methods there- 
fore are required for the determination of at least 60 chemical elements in a 
variety of uranium-base materials. Established spectrographic methods of 
analysis are unsuccessful because of interference by the complex uranium spec- 
trum in which more than 20,000 lines have been observed. 

This difficulty was overcome by separating the impurities from the uranium 
in the carrier-distillation method by converting the uranium sample to a re- 
fractory compound having low volatility (the black oxide of uranium U;Og) and 
distilling the impurities from this compound in a d-c electric arc. 

In order to sweep out the minute quantities of impurity vapors from the 
sample without volatilizing the uranium, a small amount of a volatile material, 
termed a “‘carrier,” is added to the sample. Gallium oxide was found most 
useful as a carrier and is added at a concentration of two per cent in the uranium 
oxide. When the mixture is heated by a d-c arc in a carbon electrode of special 
design, the carrier material and impurities are volatilized into the arc. The 
uranium, remaining as a residue, can be recovered readily from the electrode— 
an important consideration, particularly with some active forms of uranium. 
The light of the arc, examined with a spectrograph, provides a spectrum con- 
sisting of the simple spectrum of gallium plus the spectral lines characteristic 
of volatile impurities in the uranium. 

For quantitative determinations, carefully prepared standards of known 
composition are submitted to the same treatment as the samples. Amounts 
of impurities are then estimated by photometric measurement, or by visual 
comparison of the spectrograms of samples and standards. Determinations 
are made to an accuracy of plus or minus 10 per cent of the amount of the ele- 
ment present. For example, boron and cadmium can be observed easily at a 
concentration of 1 part in 10 million, and the concentration can be determined 
to within one tenth of this amount. The greatest sensitivity, observed for the 
detection of silver, was 1 part in 20 million. 

K. iO. 


Stand Built to Test Governors Under Service Conditions. (Compressed 
Air Magazine, Vol. 52, No. 6.)—Governors for diesel engines must be built with 
accuracy to perform the important function for which they are designed. 
These units control the amount of fuel injected into the engine cylinders by 
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the fuel pumps and come into action as soon as the engine overspeeds or under. 
speeds. In many cases where diesels drive equipment such as electric gener. 
ators, governors must be extremely sensitive to slight changes in engine speed, 
and it is for this reason that they have to undergo rigid testing and minute 
adjustment under simulated running conditions before they leave the factory, 
This applies not to one unit picked at random from a lot but to each one, and 
has led the Woodward Governor Company to design and construct a test stand 
that does this work not only in routine fashion but also with precision. 

The stand consists of a cast-iron surface plate and pedestal and is equipped 
with an air motor of the Multi-Vane Type that has a speed range from 15( 
r.p.m. to 4200 r.p.m. This unit drives the governor which, in turn, controls the 
speed of the motor by throttling the air supply. Air at from 80 to 100 p.s.i. is 
used, and the pressure must be kept constant close to the stand for satisfactory 
testing. 

Different types of governors can be coupled to the motor by means of a 
mounting plate and drive-shaft adapters, of which three are standard equip- 
ment. Depending upon whether it is keyed or serrated, the adapter is either 
attached to the governor drive shaft or dropped into a flywheel recess in the 
table top. The shaft is then inserted in the recess until the governor rests 
squarely upon the flat surface plate, something that is necessary for smooth 
running. It may also involve moving the motor up or down in order to get it 
in proper position to drive the shaft. 

After the governor has been filled with oil to a predetermined level and been 
linked with an adjustable fulcrum on the stand, the air motor is started. This 
is done by flexing the cable that connects the fulcrum and the air-control lever 
until the governor takes over. Testing involves dispelling all air bubbles by 
allowing the governor to surge for about five minutes, warming the oil by letting 
the unit run from 15 to 20 minutes, adjusting the compensating needle valve, 
and checking the response of the governor under its specified high and low 
speeds and under other service conditions. Excessive heating, vibration, oil 
leakage, and defective parts can be detected by the stand operator, and differ- 
ent measuring and indicating instruments enable him to make the corrections 
and adjustments necessary to insure a governor that will do the exacting work 


for which it is designed. 
R. H. O. 


